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ABSTRACT 


This work addresses the design and analysis of the 
Pitch, Yaw and Roll autopilot for application to the 
Bank-To-Turn (BTT) missiles. At first, the linear 
uncoupled channels were designed and analyzed according to 
the desired requirements. Utilizing the uncoupled 
channels, the linear coupled autopilots were designed, not 
including the inertial, kinematic and aerodynamic 
cross-coupling. Then, the nonlinear CBTT autopilots were 
designed and analyzed, using the linear CBTT (coordinated 
Bank-To-Turn) models, which now have coupled with 
kinematic, inertial and aerodynamic cross-coupling. The 
minimization of the above kinematic and inertial coupling 
and their effects were completed using feedbacks of 
angle-of-attack and rate of angle-of-attack in the Pitch 


autopilot. 
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TABLES OF SYMBOLS AND ABBREVIATION 


Bank-to-Turn 


coordinated Bank-to-Turn, minimum sideslip, 
positive a, 6, 188 deg 


coordinated Roll-During-Turn 


rolling moment coefficient 


Slope of curve of rolling moment coefficient, 


e vs ß 


change in с, per degree roll control incidence, 


ME 


change in C, per degree yaw control incidence, 


бү 


pitching moment coefficient 


slope of curve of pitching moment Coefficient 
= VS a 


change in Cm Per degree pitch control 


incidence, Š, 


normal force coefficient 


slope of curve of normal force coefficient CN 
ув a 


change in CN per degree pitch control 


incidence, бр 


yawing moment coefficient 


slope of curve of yawing moment coefficient, C 


n 
VS 8 


pa 
- 





pitch angular acceleration about Yg 
constant or equilibrium pitch anqular rate 
yaw angular rate about 2, 


yaw angular rate command (coordination command) 


yaw angular acceleration about Z, 


reference area for coefficients = rT ft 


Skid-to-Turn, roll attitude stabilized 


velocity component in %, direction 


velocity component in Үн direction. assumed to 
be constant 


constant missile flight path velocity 


missile velocity vector 


velocity component in 25 direction 


body-fixed roll axis, along axis of symmetry, 
positive forward 


body-fixed pitch axis, positive starboard 


vehicle axis in local horizontal direction, 
approximated as inertial axis 


body-fixed yaw axis, forms right handed 
orthogonal system with Xp and Y, 


vehicle axis in downward direction along local 
gravity vector, approximated as inertial axis 


achieved normal acceleration in 2 direction 


commanded normal acceleration in Zp direction 
achieved normal acceleration in Үр direction 


achieved normal acceleration in Ze direction 
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MY 
77 


XX 


БЕР 
LBTT 


change in Са per degree yaw control incidence, 


бү 


сһапде іп Са per degree yaw control incidence, 


OR 


side force coefficient 
slope of curve of side force coefficient Су vs 
В 


change in C, Per degree yaw control incidence, 
6 
I 


change in Cy Per degree roll control incidence, 


5n 


reference lenqth for cefficients = 2 ft. 
moment of inertia about Ya axis 
moment of inertia about Z, axis 


moment of inertia about X, axis 
autopilot pitch acceleration error qain 
CBTT autopilot coordination branch qain 


limited Bank-to-Turn may or may not be 
coordinated, positive and negative a, % 5 90 
or 45 degrees 


roll rate about Xa 


roll acceleration about Ха 


constant or equilbrium roll angular rate 


preferred orientation control 


dynamic pressure 


Pitch rate about Үр 
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achieved normal acceleration in Y,, direction 


У 


normal acceleration command from guidance 
computer in Zy direction plus anti-gravity bias 
command 

roll attitude command from guidance computer, 


zero degrees in Za direction and 90 degrees in 


Yy direction 


roll attitude, zero degrees in -Zy direction 
and 98 degrees in Yy direction 


Elevation Euler Angle, second rotation, 
(q cos$ - r sine) dt 


pitch control incidence (positive tail 
incidence produces negative pitching moment) * 


commanded pitch control incidence, бр 

yaw control incidence (positive tail incidence 
produces negative yawing moments)* 

commanded yaw control incidence, бу 

roll conntrol incidence (positive tail 
incidence produces positive rolling moment) 


commanded roll control incidence, бр 


constant or equilibrium angle-of-attack 
angle-of-attack 
ange-of-attack rate 


modified form of estimated angle-of-attack for 
autopilot coordination command 


angle of sideslip 


sideslip angular rate 
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I. INTRODUCTION 


In recent years, the application of Bank-to-turn 
guidance to missiles has generated considerable interest. 
This interest has been motivated by the drag reduction, 
which the Bank-to-turn (BTT) controlled missiles offer, 
over the conventional cruciform roll stabilized, 
Skid-to-turn controlled missiles. 

Bank-to-turn steering may provide improved performance 
for missile systems. In addition, the ramjet engine inlet 
configuration is believed to provide a greater range 
capability than other inlet configurations.  BTT control is 
required to satisfy the sideslip constraints imposed by 
chin inlets. Despite these facts, there are unanswered 
questions concerning BTT systems stability during homing, 
guidance performance, autopilot and guidance logic design 
and subsystem requirements which must be investigated 
before Bank-to-turn steering can be considered a viable 
method of control for high performance missiles. 

Many missile programs were initiated during the past 
decade to improve their capability via Bank-to-turn 
control. The results have advanced the understanding of 
the different missile subsystems. In the autopilot area, 
many types of autopilots have been found which force the 


missile to roll or bank so that the steering maneuver 
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occurs with the airframe oriented in a specified or 
preferred direction with respect to the incoming airstream. 
This entire class of autopilots may be referred to as 
preferred orientation control autopilots. 

Also, simplified guidance studies which neglect radome 
effects and assume coordinated missile motion have shown 
that CBTT can provide acceptable performance with roll 
rates that are not excessive for autopilot control. These 
Studies were made for a medium range area defense mission 
and a long range suppression mission and considered both 
high lift and moderate lift configurations [Ref. 1]. 

The present work reviews the design and analysis of 
the pitch, Yaw and Roll autopilot for application to the 
Bank-to-turn (BTT) missiles developed in [Ref, 2]. At 
first, the linear model of the uncoupled autopilot channels 
will be designed and analysed in the critical areas of 
concern regarding the response of CBTT control for the 
circular and elliptical airframes. Then, the designed 
autopilot channels will be coupled using suitable 
aerodynamic models and a measure of sideslip control is 
obtained by roll command to the linearized CBTT autopilot. 
The relative stability of the autopilot branches and means 
for improving stability will be discussed. In a later 
chapter the non-linear three dimensional model will be 


designed, using a three dimensional aerodynamic model and 
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the same control laws described in the linear studies, 
except for a minor modification to the coordinating branch 
dependence on Angle-of-Attack and also the inclusion of 
Anti-Gravity Bias. Next, the kinematic and inertial 
coupling effects between pitch and yaw dynamics will be 
studied. Transients in maneuver plane acceleration are 
caused in the above cross-coupling effects and are in the 
form of overshoots and undershoots. Coupling transients 
will be reduced by increasing pitch stability. The 
technique of feedbacks of Angle-of-Attack and rate of 
Angle-of-Attack for the pitch autopilot, will be 
considered, in order to decrease kinematic and inertial 
cross-coupling effects. Finally conclusions and 


recommendations for future study will be stated. 
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II. LINEAR DESIGN AND ANALYSIS OF UNCOUPLED AUTOPILOTS 


FOR CIRCULAR AND ELLIPTICAL AIRFRAME 


А. GENERAL 

The linear design and analysis started with uncnupled 
autopilot channels, i.e. aerodynamic, kinematic and control 
law cross-coupling between Pitch, Yaw and Roll channels 
were removed. The autopilot design technique was 
classical, using a combination of frequency respnnse and 
root locus techniques [Ref. 2], providing the range of 
required missile body angular rates, control motions and 
minimizing the influence of aerodynamic variations on 
desired responses. To achieve the desired maneuver plane 
acceleration response for the CBTT autopilot, both Pitch 
and Roll uncoupled autopilots were designed to have a time 
constant of 0.5 sec in order to achieve the desired 
maneuver plane speed of response. The Yaw channel which 
follows the Roll channel motion to produce desired 
coordination (i.e. minimization of sideslip angles), was 
designed to have a more rapid response than the Roll 
channel, having 8.39 sec for circular and 0.36 sec for 
elliptical time constant in the maneuver plane acceleration 
response. The whole design was done for the constant 
flight condition of 60 KFT altitude and mach number 3.95. 


The design discussed here was developed by Arrow [Ref. 2]. 
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In this chapter the airframe configurations, the 
aerodynamic models and control laws for the Pitch, Yaw and 
Roll channel are reviewed. The transfer functions for both 
circular and elliptical configuration are derived. An 
analysis for every channel in regard to the acceleration 
responses, the body angular rates and control surface 
deflection is made. Also, a comparison of airframe 
responses for every channel is made, considering the 
relative stability and required control surfaces. 

A general block diagram, of the channels and the 
command, which are applied to the control laws is shown in 
Figure 2.1. A number of the Figures shown in this report 


have been taken from [Ref. 2]. 


В. AIRFRAME CONFIGURATIONS 

The U.S. National Aeronautics and Space Administration 
maintains a continuing research effort in missile related 
aerodynamics. In this effort monoplanar missile concepts 
have been considered, in part, because of the relatively 
low geometric profile that can be achieved for convenient 
carriage [Ref. 1]. The two Airframe Configurations studied 
in this research were taken from (Ref. 1] and are shown in 
Figure 2.2 and Figure 2.3. The circular cross sectional 


body has a closure ratio A of 0.69 and the A... 


base/ ^max 
occurs at 68.0$ body. The Elliptical has an exact 3:1 
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elliptical cross section, and its cross section area is the 
same as the area of circular body. The span of both 
airframes is the same. The wing area and span for the 
circular concept were chosen as typical of current 
maneuvering misSiles. The wing for the elliptical 
configuration was determined by projecting the elliptical 
body on the circular bodywing planform. The resultant 
exposed wing planform then became the wing for the 
elliptical body. The tails for each design had a dihedral 
of 39° degrees as shown in Figure 2.2. In order for tail 
deflection to be compatible with complex surfaces of the 
afterbody of the elliptical configuration, the tail hinge 
line was skewed such that a 10 degrees deflection measured 
at the body-tail juncture had a resultant 7.04 degrees 
surface deflection, as shown in Figure 2.3 and (Ref. 1]. 
Comparison of the elliptical cross section concept with the 
circular cross section concept indicates the following: 

1. About 25 percent more normal force, that is nearly 
independent of angle of attack can be achieved at 
supersonic speeds. 

2. Values of the longitudinal stability parameter Cn 
are more positive, with more pronounced nonlinearities in ° 
pitching moment occuring at subsonic speeds. 

3. Levels of directional stability are increased and 


are more compatible with levels of longitudinal stability. 
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4. Noticeably more yaw control is available, although 
Suitable locations for tails on the body are more limited 


because of the geometry of the elliptical body. 


C, UNCOUPLED LINEAR PITCH CHANNEL 
l. Aerodynamic Model 

The initial phase in the design of the CBTT 
autopilots involves the design of the uncoupled pitch 
channel. The linear pitch channel dynamic model which is 
used, is shown in Figure 2.5. 

This figure is taken from [Ref. 2]. The above 
model is linearized for stability studies. The following 
three assumptions were made: 


a. Plane Хв-2 is the maneuver plane, shown in 


В 
figure 2.4 

b. Missile is trimmed in pitch, that is My=9, at 
fixed values of а, а апа $5" 

C. Missile roll rate is constant. Using the 
foregoing assumptions, and the approximation that the 
variation in forward velocity (u), is equal to zero, a 
simplified set of longitudinal equations of motion are 
constructed applicable to short period motion. These 
equations describe the two degree of freedom, short period 
mode longitudinal motion. 


Linearized aerodynamic derivatives for M-3.95 


are provided in Table I. 
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TABLE І 


Linearized Aerodynamic Derivatives (M = 
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• 0137 
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.019 


‚ 045 
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.014 


. 032 
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2. Pitch Control Law 

The pitch control laws for the circular and 
elliptical airframes are shown in Figure 2.6 as given in 
the [Ref. 2].  Lag-Leads were used to prevent guidance 
noise saturation problems. An integrator was used in the 
acceleration error branch of the circular control law to 
Satisfy the guidance requirement of zero steady state 
error. In the elliptical control law, an integrator was 
placed in the rate error path (i.e path after the adder 
(qty), as shown in Figure 2.6) for reduction of the 
acceleration response overshoot below 190 percent and a gain 
K=1.1053 was placed in series with the acceleration command 
for zero steady state acceleration error. 


A normal acceleration command n, equal to 1 дее, 


2 
С 


is applied to the pitch control law which uses measurements 
of missile body pitch angular rate (q) and pitch normal 
acceleration (n,) to determine the required actuator 


command (6$, ). The actuator is modeled as a first order lag 


P 
С 
at 188.4 rad/sec and is shown in Figure 2.5. 


3. Transfer Functions of Aerodynamic Models for 
Circular Airframe 
The equations which are represented by the block 
diagram of the aerodynamic model, of the uncoupled Pitch 


channel (Figure 2.5) are: 
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q = q Sd Са + 91.5 Са 9P | ү Iyy {vad /sec | 


Ma =~ 4S [5-3 Cusp SPt 9 ] /w | gees] (II. 
a=513[9+322n2/V] [aeg] ђе 
Utilizing Laplace Transformation, the above 


equation becomes: 


3 Gm «A + 553 Cms. d 
99: И En P) RC. 3-4) 


Iyy 
Yı= -35 Ga Cusp + Cy a) JW (II.C.3-2) 
ѕа = 513 (9+ 52215 /У) (II.C.3-5) 


Introducing equation II.C.3-2 to II.C.3-5 the last 


equation yields. 


q Sd (Са; а + 5%5 См, -Sp ) 


5 9= В EOD 
lyy 
= 322 [. 95 (513 Qu, Óp - nara ) 
pCa oo x W 
Gre?) 


Rearranging above equations, in terms of q and a, 


we obtain: 


NAT M E dics 
Sq E E - M eg)i (II.C. 3-8) 


KOS Ww. _ 515 ка 5 
-5+34+(5 + Lita) a = (- 992 0, Je ae so) 


Applying the Cramer's rule in the system of 
equation II.C.3-8 and II.C.3-9 the transfer functions q/ $5 


and a/ 6, are obtained. 
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d 
55454 („5 +[ 53 к(а5)4 (ба бб О) 
q YY T Welyy 
бр _ а 
W Iyy 
(ІІ.С.3-10) 
Апа: 
L SIKIS Case |. +Í 28329 9: 5'4 (mg 
Iyy 
Е 
B а 
С 575 35 4 Си 
52 + K. q: S Na rn o> SO Suk 
W Тү, 
езт 
Where: 
K s (C а 
М 2 
5 = 3.14159 (ft^) 
d ESO) 
q = 1659 (lbs/ft‘) 
_ во 
Е У = 804 (slug-ft ) 
W = 2525 (lbs) 
Cm = -9.06 
а 
m, = -0.08 
P 2 
EN = 9.15 
a 
€ = 0.04 
N 
Sp 
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а, = ø 
As referred to in Appendix A and Table I. After 
substituting the aforementioned values into equation 


II.C.3-18, it becomes: 


9 = 1.3361-5 — 0.459 
o (1/56 _ 
др 22.545 X 40* 57 -F 5.5633 X105 +4 ( )  (11.C.3-12) 


Introducing equation II.C.3-12 into equation 


II.C.3-2 the last equation yields: 


513 K(35) Cass a \,_ 322529 (85) 4o, Ga 


М» _ м? Iyy ` W _ 5139 5 
E = e 
Р 2 ( к. 9.5 Сне ) 57.3 q S d (ма К 
Е _ 
w Іуу 


(LI2C, 3-13) 
Substituting the above mentioned values into 


equation II.C.3-13 and minor manipulating, it yields: 


RS 
Me _ _ 406.044Х40:5* — 19.32 (gees/rad) 
- РЕТ Пе = 
др 22.545Х40 5. 534 5.3635Х40 5 +1 (II.C.3-14) 


The equations II.C.3-12 and II.C.3-14 are the 
pitch aerodynamic transfer functions for the pitch angular 
rate q about the Ув and the achieved maneuver acceleration 


in the Zp direction (n,). 


4. Design Approach and Analysis For Circular Airframe 


A Fixed Flight condition was selected for these 
preliminary performance studies of circular and elliptical 
airframes. The selected flight condition, at 68KFT 
altitude and mach 3.95, provides a sufficiently low dynamic 
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pressure, so that missile maneuvers will result in large 
enough angles-of-attack, to exercise side slip control. 
The above fixed condition was selected, in order to reduce 
the complexities of design of the CBTT autopilots. 

The uncoupled autopilot design technique was 
classical, uSing a combination of frequency response and 
root locus technique [Ref. 2] to achieve practical 
bandwidths (i.e., sufficient high frequency attenuation) 
and therefore provides the range of required missile body 
angular rates and control motions. The analysis of tne 
CBTT autopilot channels is based, on the time responses of 
manuever plane accelerartion, body angular rates and tail 
incidence angles. 

For purposes of analysis, a CSMP program was 
written (Appendix D) using the equations which are 
represented by the block diagrams of aerodynamic model and 
linear uncoupled pitch control law, as are shown in figures 
2.4 and 2.6. 

The equations of the uncoupled linear pitch 
aerodynamic model and control law are: 


a. Equations of Control Law 


(1) Acceleration Filter Equation 
X2 —————— Ma 
m td (II.C.4-1) 





Rearranging above equation II.C.4-1 one obtains: 
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Х = - '50 Х + 450 72 МЕСІ 1-2) 
(2) Acceleration Compensator Equation 


5 
- 0.381. DAL 
Y: ——————— .("-X) (1.6.4-3) 


+1) 


Utilizing inverse Laplace Transformation and 





rearranging the above equation, it becomes: 

Y +5 Y = -4.3531 ( %2,-X) -193.5X40- ("ac -X ) (II.C.4-4) 
И от state-space representation of a system, 

in which the forcing function involves derivative terms 


Meet. 5: pp. 675-673] one obtains: 


Y => ~ 4.3534 (ne, - -x) (II.C.4-5) 
Y, =-5% +6.532 - (Мә, -Х) e CEN 
ү= СЕЛЕН, 


(3) Rate Compensator Equation 





-15.6 (<< БЕ +4 ) 
(II.C.4-8) 
9р. = (XX) 
> +1) 
ы 


Using Inverse Laplace Transformation, it 


ylelds: * 


др = -0.443-Ур. - 340,58X10 (4-Х: )-2.2308 - (м - - ч) 
(II.C.4-8a) 
b. Actuator Equation (Figure 2.5): 
др = _— ‚ОР. (11.С. 4-9) 


Ja +4 
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E 
ә 


др = – 188.4 др + 183.4 др, (II.C.4-10) 
C. Equations of Aerodynamic Model 
In order to achieve the plots of the time 
response of q, noe ° in the linear uncoupled pitch 
channel, the transfer functions which are represented by 


the aerodynamic model and derived in section II.C.3, were 


used: 


О. 153 (II.C.4-11) 
ӛр 25.545 X40 $74 3.3633X40 S +14 


or 
°з э -3 B 
q (255u5X40-5^ 4.3.3633 X40-S--4) z (- 133€1-5 - 0.159). àp 
(ТС. А-а) 
Rearranging and minor manipulating equation 


ТІ.С.4-11а ІС уіе145: р 
4+449.482%10'9 + 44.5556 9 =759.2635 Ӛр -4.062-46 0 (11624212) 


Utilizing state representation of a system in 
which the forcing function involves derivative terms, as in 
section II.C.4, one obtains: 

E = X2 - 99.2635: бр (11.C.4-13) 
X = -4y.3556 X4 - 449.182X10 X2+J.789 Óp (II.C.4-14) 
а = X, (II.C.4-15) 

Also, rearranging and manipulating in the same 


way, equation II.C.3-12 yields: 


2, = > - 101901 x10" dp (II.C.4-16) 
. Ж -3 
2) =-44.3356 24 - 149.492X40 25 -1.0433X40 dp (11.0.4-17) 
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Nz = - (2, + 4.405 др) (11.C. 4-18) 
Utilizing state-space representation, all the 
aforementioned equations can be modeled in the following 


ninth order system: 


X, о | 1 О о О О О O -59.26| X, O 

К, -4.34 -0135 о оо о о о 1185 |, [о 

111 о о о 10 о о о -0||A||O l 
2, О о шн -050 о о о 045! 5) [O 

х |= о о 150 025 о O о 445 хоч, 
Y, бо о oy 1255 RS А> 


| [о о о оо 653 5 о о || с 
OR} (1231 ози O 00 -04ef -0341 -0443 -20.2 | See) baa 
llo o o oo O O а-а | РІ |0 

(Па 419) 
where the state variables are: 


X An! 


2: pitch angular rate state variables 


1! 
(section II.C.4) 
1' Yo: Achieved Body Acceleration state 


variables (section II.C.4) 


X : output of Acceleration filter 
°p 

С : input command in actuator network 
Sp : tail incidence angle 


The state X, X; Y Š of control law are shown in 


1” P 
С 
figure 2.6. 
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Executing the program of Appendix D, using an 
input step function which represents the "lGee Command” and 
trim angle of attack a, =9, the time response plots of 


pitch normal acceleration n pitch angular rate q and 


z” 
Meech tail incidence бы for circular airframe are obtained. 

Figure 2.7 shows the pitch acceleration response 
of the circular airframe due to a one Gee acceleration 
command, when the Missile Aerodynamic Models are linearized 
about zero angle of attack (a). The response has a 9.5 
seconds time constant, steady state error 9.985 and 2.5% 
overshoot, which are according to the requirements referred 
to in appendix C. 

Figures 2.8 and 2.9 show the required body 
angular rate and control surface deflection, to achieve the 
acceleration response. These responses matched those 
presented in [Ref 2]. 

5. Design Approach and Analysis For Elliptical 

Airframe 

The design technique used for the elliptical 
airframe, is the same as that of the circular airframe. 
The analysis of the Pitch channel for the elliptical 
airframe was based again, on the time response plots. 

For the above purposes, a CSMP program was 
written, as in the case of circular airframe, referred to 


in Appendix D. 
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Fig. 2.7 Pitch Normal Acceleration (n, ) vs Time (t) 
Uncoupled Pitch Channel 


Circular airframe (1Gee Command, ax =O). 
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Fig 2.8 Pitch Angular Rate (qd) vs Time (t) 
Uncoupled Pitch Channel 
Circular airframe (l1Gee Command, а.=0) 
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Fig. 2.9 Pitch Tail Incidence (5) vs Time (t) 


Uncoupled Pitch Channel 
Circular airframe (lGee Command, а, =0). 
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The equations of the uncoupled linear Pitcn 
channel for elliptical airframe, which are used in the 
above CSMP program, are: 

а. Equation of Control Law 

The control law of the uncoupled linear Pitch 
channel for the elliptical airframe is very similar to the 
control law for the circular airframe, with only 
differences in the gain of acceleration and rate 
compensators. Elliptical control law is shown in figure 
2.6., taken from [Ref. 2]. 

MI Eee rat once tE qua tion... The 
Acceleration filter equation is the same as in section 
вт. С. 4 

(2) Acceleration Compensator Equation 

Y= (-0.008/& +): (140531,,.x) (rr.c.5-1) 
Utilizing, inverse Laplace Transformation 
and rearranging II.C.5-1 it becomes: 

Y=-6Y- 530.54 X407 na, "0.9. X (11216552) 

(3) Compensator Equation 

др. _ -3.0f- (+1) 
ES - ЕНЕГЕ 5 
С-З) 
Using inverse Laplace Transformation, it 
yields: 
бр, = -383.15X10 (Y-Xj) -3.07 (Y-%) (II.C. 5-4) 
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Б. Actüator Equation (Figure 2.5). 


Кип section II.C.4,b 


C. Equation of Uncoupled Linear Pitch 


Aerodynamic Model 
The equation for pitch angular rate was 


derived, as in section II.C.3, the only differences are in 


the constants W-2475(lbs), 
20.9015, C ==-Q),.955, C =0.02 апа 
m n 
а 6 6 
р р 
С =0.18, аѕ referred in the Appendix A and Table I. 
a 


Iyy=799 (Slug-£t%), C, 


Manipulating according to the inverse Laplace 


Transformation rules and rearranging, one abtains. 


E = (3.66-5+0.638) / (-88.64X10"S7+45.95X40 +1) a 


Using state-space representation of a system, 
in which forcing function involves derivative terms [Ref. 
о рр. 675-678.] It yields: 

X4 = Х> – 44,23 др 


X; 2 11.2845X4 419.04 X40 X -15.494-dp. (11.C. 5-8) 
y RN таеке 8а) 


ев) 


Also, Manipulating with same way the equation: 


Ма z -212.441 x407 S+ 81.6 (II.C.5-9) 


др -886цХ40: 524 15.05Х40: 5 + 1 


Yielding 


-5 
Z, = 2, +528.3Х10. др (II.C. 5-19) 
2, = 11.2845-24 + 479.94 х40%, - 887.35 Op (11.С.5-11) 
Ma = — (2 + 2.936 -dp ) (11.C.5-12) 
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Utilizing state-space representation, all of 
the above equation can be modeled in the following 


eighth-order system: 





у. EE Ооо о | || о 

X. | [Hago o о Oo oo -1549||Х»| |0 

2, Оо © | о об 05131184110 

23 = | О O 4.3 013 0 0 E 8431 122 } 0 Mac 

X O 0-6о о -150 0 O-N665| X | |O 

V о о о 0 048 -60 о ||У | |-о53 

бр, 201 0394 о 090-0484 -0468 O 15.84] |Фр, | |-0.28 

бр оо о о о о 814-14 |5р| о 
[мг э БЫ) 


where the state variables are: 


X], Х-: pitch angular rate state variables 
(section II.C.4) 
Yi, Y>: Achieved Body Acceleration State 
variables (section II.C.4) 
X : output of Acceleration filter 
òp : : 
С : input command in actuator network 
бр : tail incidence angle 
The state X, X,, Yj, Sp of control law are shown in 


E 
Figure 2.6 
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Utilizing, the above state variables system, a 
CSMP program was written (Appendix D) to achieve the time 
. As input a step function 


response plots of N q and $ 


z' p 
is used, which represents a one gee command. The trim 
angle of attack (a4) is equal to zero. 

Figure 2.10 shows the pitch acceleration 
response of the elliptical airframe due to the command. The 
response has a time constant of about 9.5, steady state 
equal to zero, and overshoot o$, keeping the requirements 
referred to in Appendix C. 

Figures 2.11and 2.12 show the required body 
angular rate and control surface deflection, to achieve the 
acceleration response. 

6. Comparison of Airframe Response 

Figures 2.7 and 2.18 show that the Pitch 
acceleration response of the circular and elliptical 
airframe are approximately the same having a 0.5 second 
time constant and negligible overshoots. 

Figures 2.8,2.9,2.11 and 2.12 show that to achieve 
the acceleration response the elliptical airframe requires 
less body angular rate and control surface deflection. The 
elliptical airframe requires less control surface 
incidences, because the transfer functions have a higher DC 


gain and poles which are located at a lower frequency due 


mainly to a more neutrally stable airframe. 


cid 
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Fig. 2.19 Pitch Normal Acceleration (n,) vs Time (t) 


Uncoupled Pitch Channel 
Elliptical airframe (lGee Command, as 70) 
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Fig. 2.11 Pitch angular rate (q) vs Time (t) 
Uncoupled Pitch Channel 
Elliptical airframe (1Gee Command, а,=8) 
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2.12 Pitch tail incidence (Š ) vs Time (t) 
Uncoupled Pitch Channel 
Elliptical airframe (1Gee Command, а.=0) 
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а ш а 


р р 
proportional to pitch control moment C, and inversely 
6 
р 
proportional to magnitude of stability margin in Pitch 


The DC gain [(qS/w) .(С, RE 2 Са - CN Mares ely 
Š 


Ge., Са /с, or distance from center of pressure to 
a a 
center of gravity). Thus the nearly neutrally stable 


elliptical airframe is expected to have a higher gain than 
the stable circular airframe even though its control moment 
coefficient is slightly smaller. Higher DC gain will 
require less control surface incidence. The zeros of n2/ 85 
are directly proportional to С, and the location of the 


& 
ratio C, / Cy or the distance from the point of 
6 6 


р р 
action of tail forces to the center of gravity. 


Hence, for the elliptical airframe which has 


larger Cy the zeros of n,/ Š, 


frequency. Lower control surface incidences and higher 


are located at a higher 


frequency n,/ $, zeros of the elliptical airframe will 
simplify the design of a rapidly responding Pitch 


autopilot. 


D. UNCOUPLED LINEAR YAW CHANNEL 

The purpose of the Yaw channel autopilot of a CBTT 
autopilot is to minimize sideslip angle or provide 
Coordinated motion between Roll and Yaw channels. This was 


accomplished in two ways. First, the uncoupled Yaw channel 
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autopilot, (i.e, Roll and Pitch dynamic effect neglected) 
was designed as a regulator (i.e no guidance command and 
with rate and acceleration feedback) to help minimize 
sideslip angle. Second, to aid in sideslip control, the 
regulator was commanded from the Roll channel as explained 
in chapter III. 

1. Aerodynamic Model Control Law 

A block diagram of the uncoupled Yaw channel is 
shown in Figure 2.13, taken from [Ref. 2]. In this diagram 
the aerodynamic model and the Yaw control law are involved. 
Also, the Yaw control law without the accelerator is shown 
in Figure 2.14. These figures were abstracted from 
(Ref. 2]. 

The normal acceleration ny is not used to 
command the CBTT autopilot. Instead, it is used for the 
design and analysis of the uncoupled channel. The command 
used by the coupled system is shown in dashed lines and is 
a yaw angular rate command, ro. 

The yaw control law is governed by missile body 
yaw angular rate (r) and yaw normal acceleration п. . The 


у 
yaw control law determines the required command 6 


E 
to an actuator, which is approximated as a first order lag 
at 30 Hz (188.4 rad/sec). 


The aerodynamics, linearized about a trim а. 


angle-of-attack equal to zero. 
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ye Transfer Functions of Aerodynamic Model 


The aerodynamic characteristic of the Yaw channel 
are modeled by a second-order system, by the following 
equation, where the states are г, yaw angle and ë , 


sideslip angle: 





“= (аз4 На») - ( (ng 8 + 573 (зү: ду) (vod [sec) Dd t) 
é - 51.3 (-т + ey) (deg) 2—2) 


and also, by the equation: 

му = (9s/w)- (Cy, + ST.3 Gy * ду ) (gees) Eee 
Substituting equations II.D.2-3 into II.D.2-2, the 

last equation becomes: 

6 BET: CO (E (6:8 5x3: Cg, -5Y) ) ] (II.D.2-4) 
Taking the Laplace Transform of the equations, one 


obtains 


EE 354 . (e. 8 4 513 Gy, «dy ) mies 
hz 


45 
85 = 51.3 Ec (== (4-6 + 573 (,.4У))) 
(II.D.2-6) 
Rearranging above equations, in terms of q anda, 


they become: 


554 = HRS on.) 


Kas р - 533К35 _ 
ЗҮ + (s- R9, 6 LOT Gs OY 
ЕЮ. 2-8) 





Applying the Cramer's rule in the system of the 


equations II.D.2-7 and II.D.2-8 the transfer functions г/ 8, 














апа 8/6, аге er (44)? 
AUS IS („усу FERIEN (es Gt ©) 
ы и "ён )-5 4 UT ( ndy * Y8 УНА Суду 
by o 2 / к45 513 35d 
Ss = Су, 1.5 + Zr «c 
( W e) ( 122 ng (II.D.2-9) 
d r = 
P 53ка5 ( Y.  328»29-95d Gy 
2 —— N IS бу Ivy 
dp а 56 3954 См 
A í XE Vor )s- = 
ҮУ (II.D. 2-19) 
(Б7з) (5325) 
EE ne - = 0.48 
V 
q - 1650 (lbs/ft?) 
5 = 3.14159 (£t* 
122 = 819 (slug- ft“) 
м = 2525 (lbs) 
С = -8.025 
B 
En 
5, = -0.050 
C 
Ув = -0.065 
EY 
$ = 
" 0.021 
a = po 


As referred to in Appendix A and Table I. 
Substituting the aforementioned values into equation 


(II.D.2-9), it becomes: 60 





v .  2.005.$ + 0.46 
бр,  -5u.54X10-S-2.8X10 3 5 t 
Introducing equation II.D.2-10 into II.D.2-3 and 


(ilsec) (TI: D.2-11) 


rearranging the last one, the transfer function of yaw 
normal acceleration in terms of the yaw control incidence, 


is obtained: 





513 x (95)? қ 323229 (84) d 
— go JS HERTZ E ay E 
My - | w у “Ye Iyy W now E 5385 Gy 
Sy «141 чи) 5- 513 5 Sd Gae ы 
м Туу 


№. 2—2) 
After substituting the values into the equation 
(II.D.2-12), one obtains: 


му _ -43.44Н4 SY 47.16 (es rad 
dy -5ы.5ҢХ10:55-2.8Х10:5 +1 


(II. D.2-13) 
The above transfer functions are used for finding 
the time responses of the yaw angular rate r and the yaw 
normal acceleration. 
Using inverse Laplace transformation, the equation 
ШЕП; 2-11 and II.D.2-13 become: 


e -3 о 
“ +54348Х10- У ~ 19:3535 7 = -36.464- ду - 2.45- Sy (о 14) 


г © 


My +51.338X40° y -183335 ny = 2.387 Sy -325.629-dy (II.D.2-15) 


The equations II.D.2-14 and II.D.2-15 involve 
derivative terms in the forcing function. Using rules of 
state-space representation of a system [Ref. 5: pp 


675-678], they turn: 
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Ха - X3- 36-161. Óy : 2 Lo) 
X, = 18.335-X, - 51,338X10-X, -363.044 X40-9Y  (11.D.2-17) 


Y = XI : unm 18) 
2, = 22 - 122.544 x10-0y E. 
2) = 48.335-24 - 51.338X40 22 - 291.851.8y  (11.D.2-20) 
Му = 24 + 2.397. ду (стар. 2=21) 


E. Equations of Yaw Control Law 


As mentioned before, the yaw control law is 
governed by missile body yaw angular rate (r) and yaw 
normal acceleration N, and it is used to determine the 


required command Š, to the actuator. The yaw control law 
at the flight SLT ER of interest (i.e, Mach 3.95, 60kft) 
is shown in Figure 2.14, taken from [Ref. 2]. The rate 
compensation determines the high frequency attenuation and 
is used to minimize aerodynamic variations on the quality 
of regulation. The acceleration compensation determines 
the acceleration bandwidth via the time constant of the 


acceleration response due to a step command of acceleration 


at ny • 
С 
а. Acceleration Compensator Equation: 


E 0.31946 т (II.D.3-1) 
NE CSS (mes y) 


Inverse transforming and rearranging the 
II.D.3-1, it turns into: 


Y = -5Y + 1.5933 ( y - ус) (II.D.3-2) 
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b. Rate Compensator Equation: 


5 
Sy, = 4.85 (35 +1) (y+r) 09858) 


Utilizing inverse Laplace Transformation 
and introducing equation II.D.2-16, the equation II.D.3-3 


turns into: 
бу. = 6.435 (У+Х4) +4. 85 (У+Х4 ) За) 
с. Actuator Equation: 


jm : 1 us (II.D.3-4) 
uma t! 
Inverse transforming and rearranging the 
equation II.D.3-4, it turns into: 
ду = - 188.4. бу + 188.4 Se (II.D.3-5) 
Using state-space representation, all the 


aerodynamic model and control laws can be modeled in the 


following seventh order system: 


Xi O | о О О О -36.16 X4 О 
Х | |18335-054 © о о о -0.863{ | X2 O 
2. 0 О O 1 О О -04226 21 О 
2, || о о 15335-0584 o o 8190“ |2 |410 | 
Y O о 15% O -5 о 3399 |Y -1.591 
$, | 14.95 0.485 045 0 2425 О -1828| 19% | 10.798 
by o о о о о 1884 -4884 
(II.D.3-6) 


where the state variables are: 


X), X^: pitch angular rate state variables 


(section II.D.2) 


63 





Yi: Ү,: yaw acceleration state variables 


(section II.D.2) 


Y : output of acceleration compensator 
Sp 
E : input command in the actuator 
бр : yaw tail incidence 
The state Ху, Y $, of control law are shown in 
€ 


Figure 2.14. 


4. Design Approach and Analysis of Circular Airframe 


The design technique used for the Yaw channel is 
the same as the technique used in the Pitch channel. The 
analysis of it, for circular airframe is based, on the time 
responses of the yaw angular rate r, the yaw normal 
acceleration m and yaw tail incidence m 

For purposes of analysis, a CSMP program was 
written (Appendix E) using the state-variable system which 
is represented by the block diagram of the yaw aerodynamic 
model and control law. 

Executing the referred program, using as input a 
step function, which represents the "one Gee Command" and 
trim angle of attack а.=0, the time response plots are 
obtained. 

Figure 2.15 shows the yaw acceleration response of 
the circular airframe. The response has a 0.39 seconds 


time constant, 7% overshoot and steady state error 0.018 
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NY (GEES) 





Fig. 2.15 Yaw normal acceleration (ny ) vs time (t) 
Uncoupled Yaw Channel 
Elliptical airframe (lGee Command, a =0) 
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which are according to the requirements referred in 
Appendix C. 

Figures 2.16 and 2.17 show the required body yaw 
angular rate and control surface deflection, to achieve the 
acceleration response. Again these results matched those 
of [Ref. 2]. 

ЗС хегодипаті е Transfer.bunction for Elliptical 

Airframe 

The transfer functions for yaw angular rate and 
yaw normal acceleration for the elliptical airframe, were 
derived in the same way as in section II.D.2. The 


only differences are in the following constants 


М = 2475 (lbs) 
122 = 853 (slug-£t? 
C = 0.024 
ng 
^n, = -0.042 
Y 
Су, = 0.016 
ү 


As referred to in Appendix A and table I. 


a. Transfer function of yaw angular rate r. 


V _ -13435-5-0.059$8 1 
oS —— (ser) (11.D.5-1) 
у 59.858 10 .5 4 2.593 10: 5 + 1 
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Fig. 2.16 Yaw angular rate (r) vs Time (t) 
Uncoupled Yaw Channel 
Circular airframe (1Gee Command, а.=0) 
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Fig. 2.17 Yaw tail incidence (бу) vs Time (t) 


Uncoupled Yaw Channel 
Circular airframe (1Gee Command, a ,70) 
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b. Transfer function of yaw normal acceleration пу 


-3 
Мү - 414.633 Y40- $—- 1.1 gees 


" = — DS 
бу 59.938 10 3 574 2.593Х40-5+4 V Yad ) (11.D.5-2) 


Using inverse Laplace transformation, the equation 
О. 5-1 апа II.D.5-2 becomes: 
Z13333 X10: Y446.1117- Y =-29,243.ду- Q99.36X40-dy (IT.D.5-3) 
Му 343.333Х1074у 4163ШЗ-Ту -1.0453-2у-48.663-6У7 (11.0.5-4) 
Utilizing State representation of a system in 
which the forcing functions involve derivative terms, as in 


section II.D.2, the above equations yield: 


у. = Хо- 29.221 ‚дү (II.D.5-5) 
E - -ACTHT X47 03.333 40- Xo 266.3148 Y40-0y — (11.D.5-6) 

= X (II.D.5-7) 
4 - 22- 93.013X40. dv (II.D.5-8) 
2, Z IC HI}. 24 ~ ц3,333Х 10: 25 –150.6633 · ду (11.0.5-9) 
My 24+ 1.954% -dy (11.0.5-10) 


The aerodynamic model for elliptical airframe is 
shown in Figure 2.13. 
6. Yaw Control Equation for Elliptical Airframe 
The control law of uncoupled linear Yaw channel for 
elliptical airframe, is similar to the control law for 
circular airframe. There are differences in the gains of 
the acceleration and rate compensators, as is shown in 


Figure 2.14. 
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a. Acceleration Compensator Equation 
Y= (0.8393/0.25-5 +4) ' (-му, + tv) (11.D.6-1) 
Utilizing inverse laplace transformation and 
rearranging the equation II.D.6-1, it turns into: 
Y = -4 y + 3.3574 · (-ту +) (II. D. 6-2) 
b. Rate Compensator Equation 
gy - 6.08. ett) (Y-v) 
5 ео. 63) 
Inverse Laplace transforming and introducing 
the equation (II.D.5-7) into the equation (II.D.6-3), it 
Barns into: 
бу, = 0.608 (Y+X;) + 6.08 (Y+X1) (ШЕ b. 6-4) 
C.  Actuator Equation 
Same as the equation II.D.3-5. 
Using State-space representation all equations of the 
aerodynamic model and control law, can be modeled in the 


follow seventh-order system: 





X O 1 О О О o 292 |Y. o 

Xa | 492-0435 o о о о 0.241 |Х2 о 

, O o 0 T 0 о -0.85 24 О 

5 || о o -en-o о о 06|. |2 |+| 0 |, 
у о о 356 © ^H О 6.553 Y -3,36 

бу, 6.08 0.608 204 0 -2432 O 3.984 oc -2. 04) 

y о о о о о 1884 -1884| | V Л) 


(II.D.6-5) 
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where the state variables are: 


Xi, X^: yaw angular rate state variables 
(section II.D.5) 
21, 22: yaw normal Acceleration state 
variables (section II.D.5) 
Y : output of acceleration compensator network 
a 2 input command in the actuator 
Sy ‚ yaw tail incidence 
The state Xi; ү, Š y of control law are shown in Figure 


С 
2.14. 


2% Design Approach and Analysis for Ellipticai 


Airframe 

The design technique and analysis for elliptical 
airframe is the same as that, in section II.D.4. Again, 
for purposes of analysis, a CSMP program was written 
(Appendix E) using equations of aerodynamic model and 
control law. 

Figure 2.18 shows the yaw acceleration response of 
the elliptical airframe. The response has 9.35 second time 
constant, overshoot about 6$ and steady state error 9.01 
which are according to the requirerments, (Appendix C). 

Figures 2.19 and 2.20 show the required body yaw 
angular rate and control surface deflection, to achieve the 


acceleration response. 
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Yaw Acceleration (пу) vs Time (t) 
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Uncoupled Yaw channel 
Elliptical airframe (lGee Command. а,=8) 
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Yaw Angular Rate (r) vs Time (t) 
Uncoupled Yaw Channel 
Elliptical airframe (lGee Command, а =й) 
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Fig. 2.20  Yaw Tail Incidence (6,) vs Time (t) 
Uncoupled Yaw Channel 
Elliptical airframe (lGee Command, а.=0) 
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Е. LINEAR UNCOUPLED ROLL CHANNEL 

Generally, the purpose of a roll channel autopilot is 
to command the ¢ signal, in order to roll the missile 
through an angle $, measured from the vertical axis. 

The roll channel of a coordinated bank-to-turn 
autopilot is commanded to roll the missile so as to put the 
preferred maneuver direction of the missile in the 
direction of the guidance acceleration command, while the 
pitch channel acceleration is commanded to produce the 
total magnitude of the guidance acceleration command. 

The desired maneuver plane acceleration should be 
attained as rapidly as the achieved body-fixed pitch 
acceleration. To accomplish this, the uncoupled roll 
channel autopilot was designed to have the roll angle time 
constant equal to the time constant of the normal 
acceleration achieved by the uncoupled pitch channel 
autopilot. This is not necessarily the optimal 
relationship betweem the two time constants. 

l. Roll Aerodynamic Model and Control Law 

A block diagram of the uncoupled Roll channel is 
shown in Figure 2.21,taken from [Ref. 2]. In this diagram 
the aerodynamic model is involved. Also, the roll control 
law is shown in Figure 2.22. These figures were taken from 
mer. 2]. 

The roll control law is commanded by roll angle ve 


iS governed by roll angular rate p and roll angle ф апа 
75 
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determines the required command 5n to an actuator which is 
approximated as a first order lag 28 30Hz (188.4 rad/sec). 
The aerodynamics, linearized about a trim angle-of-attack 
(aa) equal to zero. 
2. Equations of Aerodynamic Model for Circular 
Airframe 


The aerodynamic model is represented by the 


following equations: 


Р = (J Sd lixx) Сер + Og (хас |5еС) (11.Е.2-1) 
Ф - F |5ес) (II. E.2-2) 
where: 
q = 1650 (lbs/ft‘) 


s = 3.14159 (ft^) 


d = 2 (ft) 
Тхх = 49 (slug-£t”) 
и“ = 0.031 

R 


As referred to in Appendix A and Table I. 
substituting the above values in equation II.E.2-l, it 
yields: 

P = 8.0346 -dp (II. E. 2-3) 
3. Equations of Aerodynamic Model for Elliptical 
Airframe 


The aerodynamic model equations are similar to 


those for the circular airframe, with the only differences 


T3 





in the constant 1, ,=110(slug-£t”) and С, 2505 023. 


5n 
Hence they are: 


p = 9.1643 ‚др (rad sed) 
Фф = Р (коо іѕес) (II.E.3-2) 


4. Equations of Control Law for Circular and 
пресса Ат гегате 


The yaw control law for circular and elliptical 


ИЗ) 


airframe is shown in Figure 2.22. 
а. Roll Angle Compensator Equation 


И = — ( 4-6) 


+1 (II.E. 4-1) 





$ 
inverse Laplace transforming the equation II.E.4-1 


it becomes: 


Х =-8Х+11.6 ( be- $) (5,12) 
b. Rate Compensator Equation 
[== 
y. 21595 (Ren («-P) ER. 


(++ 1) 
Utilizing inverse Laplace transformation, and 
minor manipulating the equation II.E.4-3, one obtains: 
Y 7 -5-Y+ 50.333X4 * (х- р) + 0.3555 (Х-Р) (ТТ.Е. 4-4) 


с.  Psuedo Differentiator Equation 


"Nie 0.01303 


5 
6 


‚р (ТТ.Е. 4-5) 
+4 
Using inverse Laplace transformation and 
rearranging the equation II.E.4-5, it becomes: 
X = -6% +0.0432.Р (ТТ.Е. 4-6) 
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а. Equation of Actuator Compensator 
This network of roll control law is different 
r circular and elliptical airframe. 
(1) Circular Airframe 
= E) CY-X4) / (+ + 4) Е.Е. 4-7) 
AL, inverse ШІ transformation and 
rearranging the equation II.E.4-7, it yields: 


Ф. = -455р + 0.4565 (У-)+15(1-Х) (11.8. 4-8) 
(2) Elliptical Airframe 


КЕ ы (£,+ 1) ey. (II.E.4-9) 


Inverse transforming and 


rearranging the above equation, it becomes: 


SR. = - 155g, + 0.569393 (Y-X, )+62.536 (Y-xı ) 


e. Actuator Equation 


ба, = — — дес (deg ) (ІІ.Е.4-11) 


488. жі 


Utilizing inverse Laplace transformation, the 


equation II.E.4-11 becomes: 


bp - -1884.óp + 10195.32 dg, (del (re. 4-12) 
Utilizing state-space representation, all the 


aforementioned equation of the aerodynamic models and 
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control laws for the circular and elliptical airframes, can 
be modelled in the state variables systems given below: 


(1) For the circular airframe 


5 О О О О О О 8.0345 р © 
ф { 0 О О О О О ф О 
Х O $6 -8 O о 0 о Kal JI 
У |=|-045 05% 0 -5 о 0353-0404 |*| Y |+|-0.889 Ф 
X % 0 O O -6 07 | |Х4 0 
бес | Loto 0421 0.483 Н.м9 -щие -5 -0.033 | [Se | |-0.12 
ӛр O © о о о 10932 -1394 | [og О 


Е) 


(2) For the elliptical airframe 


Р о о о д о о ¿ie pp О 
ф Í 0 ОР Eoo O ф О 
X о 6 о о о о ||х| [6 
ү |= [-025ь 080 О -5 0355 0 -одощ У | т |-0,856 фе 


X Do oo c O 633 |X1 | O 


бас -4.243 -0.504 0.208 3.441 59.448 -15 -0.58% Оре -0.504 
M о о O о о 1095 - |Ф | |о 
(II.E.4-14) 
where the state variables are: 
P : roll angular rate 
$ : roll angle 
X : output of roll angle compensator 


Y : output of rate compensator network 


QE 





output of pseudo differentiator 


6, : input command in the actuator 


Sa * roll tail incidence 


1! 5; are shown in Figure 
С 


The state $,X,P,Y,X 


522. 


5. Design Approach and Analysis of Uncoupled Linear 
Roll Channel for Circular and Elliptical Airframe 


The design technique and analysis is the same as 
that which was used for the Pitch and Yaw channel. 

For SS of analysis, two CSMP program were 
written (Appendix F) using the state variable system 
represented by the aerodynamic model and control law, 

Block, diagram. 

Figures 2.23 shows the roll angle response of both 
elliptical and circular airframe. The time constant of the 
roll angle response is 0.55 seconds, the overshoot 3$ and 
the steady state roll angle error equal to zero. 

Figure 2.24 shows the required body roll angular 
rate for both elliptical and circular airframes. Only the 
roll tail angular deflection figures 2.25 and 2.26 is 
different for the airframes, due to the method of 
compensating for a reduction in the elliptical aerodynamic 
roll gain. This is desirable for the aerodynamic roll gain 


(57.3qsdC, /I,,) to be as large as possible to minimize 
6 
R 
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Roll Angle (9$) vs Time (t) 
Uncoupled Roll Channel 

Circular or Elliptical airframe 
(1 RAD Command, as 79) 
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Roll Angular Rate (p) vs Time (t) 
Uncoupled Roll Channel 

Circular or Elliptical airframe 
(1 RAD Command. а =й) 
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Fig. 2.25 Roll tail incidence ($,) vs Time (t) 
Uncoupled Roll Channel 


Circular airframe (1 RAD Command, 
a 7-9) 
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Fig. 2.26 Roll tail incidence (Sp) vs Time (t) 
Uncoupled Roll Channel 
Elliptical airframe (1 RAD Command, а,=8) 


56 





control surface motion. 

large C, as possible to 
aerodynamic control cross 
has a considerably larger 
much smaller roll inertia 


E. е 
К 


It is also desirable to have as 


minimize the effects of 
coupling. The circular airframe 
aerodynamic roll gain due to a 


and a larger control derivative 
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III. LINEAR DESIGN AND ANALYSIS OF COUPLED AUTOPILOTS FOR 


CIRCULAR AND ELLIPTICAL AIRFRAMES 


A. GENERAL 

The linear design and analysis technique began with 
uncoupled autopilot channels. The uncoupled autopilot 
design technique was classical, using a combination o£ 
frequency response and root locus techniques [Ref. 2], to 
achieve practical bandwidths (i.e, sufficient high 
frequency attenuation) and in turn provides the range of 
required missile body angular rates and control motions, as 
mentioned in section II. 

In this section, CBTT control laws are analyzed to add 
control coupling for coordinated missile motion. А measure 
of sideslip control is obtained by applying a roll angle 
command to the linearized CBTT autopilot. The relative 
Stability of the autopilot branches and means for improving 
stability are discussed. An examination of the autopilot 


sensitivity to aerodynamic cross-coupling is made. 


В. AERODYNAMIC MODELS FOR THE CIRCULAR AIRFRAME 

The aerodynamic models of pitch and roll-yaw channels 
for the circular airframe are shown in Figures 3.1 and 3.2, 
taken from [Ref. 2]. The models are linearized and were 


developed for stability studies in the frequency domain. 
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The assumptions for linearization are the same, as those in 
K ion II.C.1. 

Pitch and Lateral channels are coupled via constant 
missile role rate Poe Also, the linear model has inertial, 
kinematic and aerodynamic cross couplings, which will be 
discussed further, in the non-linear model. 

l. Pitch Channel 

The linear aerodynamic model of the pitch 
channel consisted of the uncoupled model mentioned in 
chapter II, which is coupled with the Roll-Yaw channels via 
constant roll rate (Р), as shown in Figure 3.1. 

In this work, the Pe was set to zero, because 
when Ре-0, the Qe (i.e equilibrium pitch rate) has been 
found to have negligible influence in the lateral model 
(i.e roll-yaw). 

The equations, which represent the pitch 
aerodynamic model given in [Ref. 2] are the same as the 
equations of the uncoupled linear model, mentioned in 
chapter II. 

2. Lateral Channels 

The linear lateral aerodynamic model consisted 
of the uncoupled roll and yaw aerodynamic channel, coupled 
via the linearized aerodynamic derivatives Co и С, апа 

8 бу 
C. . Also, it is coupled with the pitch dynamic model via 


SR 
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the Pe, as previously mentioned. The lateral dynamic model 
is shown in Figure 3.2, as given in [Ref. 2]. 

To find the equations which represent the 
lateral dynamic model, Pe and Qe are set to zero because of 
their negligible influence in the lateral channel. 


The aforementioned equations are as follows: 


y= -S-E / үу (gees) (II.B.2-1) 
E42 6: Qa + Сур, Óy (deg) (T Beo) 
where: 
ç 2 9533935-4 | 12а ( degisec”) NS 
Ь= 6. (лв + (ме „Фр+ Сту. ду (дее) $59 9E PPP S 
P= 53.39 Sd Ey [ Ixx (ес [зес2) (III.B.2-5) 
ф- РІ5%3 (radlsoO (111.B.2-6) 
533 Ces, .ду + Cese- ÔR (deg) (III.B.2-7) 
Ey= Ea + Сед. 6 (deg) (III.B.2-8) 
4- PA[sr3- X*W-(q3 Ki] w) (deglsed (ттт.в.2-9) 
where: 
k - 0.48 
S - т(ға)2 
а = 2.0(£ft) 
q = 1650 (1b£/ft2) 
w = 2525 (lbf) 
Ixx - 40.0 (slug-ft?) 
Izz = 819.8 (slug-ft2) 
u = 0.22 


Sa 
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= 0.082 
E 
n 
Š = 
R 0.018 
С 
ПВ = -0.019 
С 
n 
5) -- 
у 0.053 
Cy 
3 Че 
ү 0.016 
С, 
В ш -0.009 
Cy 
6 = 
R 0.035 
À = 10° 


As referred to in Appendix A and Table I. 
Substituting the values into the above equations, 
they yield: 


Ф 


y 


-43.934 -6 + 13.954 др - 38.81.0У (deglsed) (IDEO NOD) 
- 133.66 -8 + 519.183-02 -234.618.4y(deglsed) (111.8.2-11) 


р 
& =0.01345 P-Y-0.0808 8+0.17-0Y (degkec) (T17-B- 2-12) 
ф Plsy.3 (тоа [ес (ІІІ.В.2-13) 


i 


С. CBTT AUTOPILOT CONTROL LAWS FOR CIRCULAR AIRFRAME 

The control laws, which were used by the CBTT 
autopilots of the circular and elliptical airframe, are 
Shown in Figures 2.6 and 3.3, as given in [Ref. 2]. 

The pitch control law for the circular airframe is the 


Same as determined in the uncoupled pitch channel study 
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section II.B.2, and therefore, the equations which 

represent it, are the same as in section II.B.3., except 

the equation for pitch tail incidence angle which is: 
др= - 433.4 др +40 +95. 3 ·др, (деда бес) (QUEUE (Cic dl 

The uncoupled roll and yaw laws are coupled via a 
cross-coupling branch, shown in solid line in Figure 3.3, 
in order to consist the linear coupled lateral control law. 

In the roll actuator of the lateral control law, the 
frequency has been changed from 118 rad/sec to 60 rad/sec, 
for improving of autopilot stability. 

The autopilot cross-coupling branch between Roll and 
Yaw channels has been added to provide coordinated motion. 
Coordinated motion or zero sideslip angle is achieved by 
directing the body fixed pitch axis of the missile of the 
missile velocity vector, so that there is no component of 
missile velocity along the body fixed yaw axis of the 
missile. Figure 3.4, as given in [Ref. 2] shows the 
attitude of the missile body with respect to its velocity 
vector V. 

When commanding an upward maneuver (i.e, Ф =0), the 
missile body moves upward with its pitch axis directed at 
the velocity vector until it reaches the desired maneuver 
level or angle-of-attack. No roll motion is required to 


maintain coordination for this maneuver. 
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Fig. 3.4 Coordinated Missile Motion for 
Coordinated bank-to-turn (CBTT) 
and limited bank-to-turn (LBTT) 
Control policies 
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For maneuvers in the %=45 or 90 degree directions, the 
missile body rolls the desired % and moves its pitch axis 
directed at the velocity vector, until it achieves the 
desired maneuver level. For maneuvers with ® higher than 
99°, the missile must roll about the velocity vector, while 
the Yaw channel directs the pitch axis towards the velocity 
vector for minimum Sideslip angle. 

The missile is forced to roll about its velocity 
vector, when the desired maneuver direction 1S in the 
negative angle-of-attack direction, because it is desired 
to avoid negative angles-of-attack. These maneuvers are 
illustrated in Figure 3.4. 

The coordinating command ra, is a yaw angular rate from 
a rate gyro. The command is equal to - Рас: where a,=a 
(Figure 3.3) in the coordinating branch is the equivalent 
angle-of-attack, which is exactly equal to angle of attack 
(a), in the linear studies. The equations which represent 
the lateral channel are as follows: 

1. Roll Channel 

a. Roll Command Compensator Equation 
Same equation as in section II.E 
X2 9x 44-6 (de-$) (III.C.1-1) 
D. Roll Rate Compensator Equation 
0.151: ( +1) 
= (х- Р/аз) 





т.е. 1—2) 
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Inverse Laplace transforming and rearranging 
the equation (III.C.1-2), it becomes: 


Y, 559 4+ 0.0503 (X-Plst3) + 0.755 (x - p/sr.3 ) Вет. 1-3) 


с. Psuedo Differentiator Equation 


ТЕ 0.015 - 5 ( p/ 57.3 ) ГІЛІ 7.1224) 


> 4 
C + 


Utilizing inverse Laplace transformation and 





minor manipulating the above equation, one has: 
Х, = -6Х, + 0.048 (Р/513) ШОО О; 
d. Actuator's Compensator Equation 
The frequency of the compensator was changed 
Ғгоп 110 rad/sec to 68 rad/sec for purposes of stability, 
as referred in section III.C. The equation represents the 


network iS: 


- + | 


Ec (Y- X0 (IFI GC =O) 


uA d 
15 


Using inverse Laplace transformation, the 
equation III.C.1-6, one obtains: 


ЈЕ 2 2122 др. + 0.25 ( Y-X,) +45 (у-Х‹) (TITI 7) 


e. Actuator Equation 


dp = 245 de, (4а) (III.C.1-8) 


E 
188.4 


Utilizing inverse Laplace transformation and 


rearranging the equation III.C.1-8, it yields. 
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ӛр - - 188.4 da + 10395.32 бес ОЙГО о) 


All the referred to above state variables, are 
shown in figure 3.3. 
2. Yaw Channel 
a. Acceleration Compensator Equation 


> 0.32 (те С 22-1) 


N 
025+ 1 М 


Inverse Laplace transforming, equation 
III.C.2-1 becomes: 
У = -5У +16 Ny (pi 252) 
b. Actuator Compensator Equation 
ду, = (445(-5--н) [$). (Ү- 0.958 а.РІ5.3+1) (I11.C.2-3) 
Utilizing inverse Laplace transformation and 
rearranging the equation III.C.2-3, it yields: 
ду, z 0.485(Y+Ylsr3-0.458-a- P/51.3)+ 4-25(Y+ rls1.3-0.458-3-P/51.3) 
Е.С. 2—4) 
с. Actuator Equation 
ду = ~ 138.4 -òy + 40795.3 - dy, ne 
Making use the state-variable representation, 
the above equations, of the aerodynamic models and control 
laws, can be modelled in the following, state variable 


equation system: 
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variables are: 

: yaw angular rate 
roll angular rate 

I roll angie 

>: Sideslip angle 

: output of roll command compensator 
output of roll rate compensator 


: output of pseudo differentiator 


: input command in roll actuator 


: roll tail incidence 
: output of yaw acceleration compensator 
: input command in yaw actuator 


: yaw tail incidence angle 


the referred before state variables are 


Shown in Figures 3.2 and 3.3. 
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D. ANALYSIS OF LINEAR CBTT AUTOPILOTS FOR CIRCULAR 

AIRFRAME 

In the analysis of the linear CBTT autopilots it is 
assumed that the missile is initially in the desired 
maneuver plane, is trimmed at ten degrees angles-of-attack, 
the equilibrium roll rate Pe and equilibrium pitch rate Qe 
were set equal to zero. 

For purposes of analysis, a CSMP program was written 
(Appendix G) using the equations of the aerodynamic models 
and control laws of the linear CBTT autopilots for the 
Specular airframe. 

Figure 3.5 shows the sideslip angle for the linearized 
CBTT autopilot, having a maximum value of 25222 Since the 
model is linear, and therefore the magnitude of the 
Sideslip is directly proportional to the magnitude of the 
input command, a one radian roll angle command was used for 
convenience. 

The coupled autopilots have a roll angle response which 
is the same, as the uncoupled roll channel response, 
differing slightly in the overshoot, as shown in Figure 
3.6. 

The aforementioned values and plots are exactly the 
same as the plots given in the [Ref. 2] therefore its 


results are verified. 
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Fig. 3.5 Sideslip angle (g) vs Time (t) 
Linear CBTT Autopilot (1 RAD Roll 
Command); Circular airframe (а -10 deg, 
Pe=Qe=9 deg) 5 
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Comparing the state variable form system with these 
previous discussed in the section II.D and II.E, this is 
the twelfth order system because it includes roll and yaw 
autopilots and the matrices A,B are different than the A,B 
matrices at the uncoupled roll and yaw autopilot, because 


of the existance of the coordination branch. 


Е. AERODYNAMIC MODELS FOR ELLIPTICAL AIRFRAME 

The dynamic models for elliptical airframe are the same 
as those of circular airframe and are shown in Figures 3.1 
ШІП 35.2. 

The differences of the circular and elliptical models 
are in the following constants and linearized aerodynamics 


derivatives: 


W = 2475 (156) 
Ixx = 118 (slug-ft*) 
122 = 799 (slug-ft*) 
С 
Y — 
C 
Ys = -0.054 
En 
5а = 0.014 
En 
в = 0.024 
С 
п — => 
5, = -9.039 


104 





“ы 010 
MEL 
C 
ба - -0.027 
E. 3 moms 
R u: © 


as referred to in Appendix A and Table I. The equation 


which represent the aerodynamic model, are as follows: 


+ = 49.0069.6 410.5213. 0g - 29.3262. y $95 48-958 
$ = -145.4104- 6+424.1094-др - 54.04. ду (III.E.2) 
$ - Pls+3 (TOTIES) 
Ё = 0.13452 p - Y +0.05428F-8-0.01508.0Y — (III.E.4) 


F. CBTT AUTOPILOTS CONTROL LAWS FOR ELLIPTICAL AIRFRAME 

The control laws, which were used by the CBTT 
autopilots of the elliptical airframe are the same as those 
of the circular airframe and are shown in Figure 2.6 and 
БЕ 3. 

The differences which exist in the control laws of the 
two airframes are the gain Кур of the coordinated branch, 
and the gains of the actuators. 

Hence, the equations that represent the control laws 
are: 

1. Pitch Channel 

The linear pitch control law equations for 
the elliptical airframe, are the equations of the Pitch 


uncoupled channel, as given in section II.C.5. 


5015 





2. Roll Channel 


Same equations as in par. III.C. 


a. X =-8х+ 44.6 (фе-Ф) (ZI.6.1-1) 
b. Y, 2 -5Y,t 0.0505 (X-PIsr.5) c 0.155 (x- P/57.3) (11.C. 1-3) 
с. Х, = -6Х, + 0.048 (Р/51.3} (11.C.3-5) 


d. Actuator Compensator 


(ТЕ... 2-1) 
te Att Ge 4) cy- Y 
uel 


Utilizing inverse Laplace transformation and 
rearranging the equation III.F.2-1, it becomes: 
др, = “155g. + 1.0425 (Yy-X ) + 62.53 (М,-У,) (ТЕШЛЕ. 2=2) 
е. бр = ~ 183.4 ӛр + 10.4955 Х16 бе, (тес. 3) 
3. Yaw Channel 


a. Acceleration Compensator Equation 


wre 
ос 
Д 


Inverse Laplace transforming and rearranging 
the equation III.F.3-1, it yields: 


У = -ЧУ + 3.556 Ny (TTI. F. 3-2) 


b. Actuator Compensator Equation 


Bye = CBS ESOS ( y+ rista -2 P/5F3) 
(ITI. E, 3-3) 
Utiliing inverse Laplace transformation and 


rearranging the above equation, one тес 


v ap 
Sy, = 0.608 Ст = + 6:08 (Y* sz — = 2 (III.F.3-4) 
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the equation which represent the 


control 


с. 


Same as in par. 


Actuator Equation 


РЕН ве се 


Utilizing the state-variable representation, 


laws, 


can be modelled 


twelfth-order system: 


in 


aerodynamic model and the 


the following 


ү ЕЕ: 2 50 о о О O 0 39929| |” 9 
P о о оша о о о оши о о ww |? 0 
Ф Шо о о о со о о о о о о ф 0 

g Он 0 00580 0 0 0 0 0 о Hori] 16 0 
X ЕНЕ 2 ПП о -8 о о о O О ООо x (1.6 
{|1 0 -0192-0986 015 035 5 о о 00 о о 004T|* Y [4 | 036 | - 
| | o 0 o -nt o о 6 о 968 о o0 4R] |X} уо 
Be O 0.0137 -0926 1199 Ox 5144 -568 -15 4024 о o Чаш | || | 0.926 
1 о оо о о о © 53-1834 о о O OR О 
Y ЕЕ) 20550 о о о о -4 о 905| |y 0 
By} | 0106 -oons o oise o о О © os 534 0 мн! | 0 | 
| | О О O 0 0 0 0 0 О O 40795.3-488-4 | 15ү 0 


(ME BE 3-5) 


where the state variables are: 


y 
~ 


р 


yaw angular rate 
roll angular rate 
roll angle 
side slip anqle 
output of roll command compensator 
output of roll rate compensator 
output or psuedo-differentiator network 


input command in roll actuator 


roll tail incidence angle 
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ER output of acceleration compensator in yaw 
control law 


т; input command in yaw actuator 
DC: yaw tail incidence angle 


All the state variables are shown in Figuxe 3.2 and 
3.3. The above system is different in the matrix A and B 
from the system discussed in chapter II for uncoupled Roll 
and Yaw autopilots, because of the coordination branch 


equation. 


G. ANALYSIS OF LINEAR CBTT AUTOPILOTS FOR ELLIPTICAL 

AIRFRAME 

In the analysis of the linear CBTT for elliptical 
airframe, the same assumptions, which were made for the 
circular airframe, are valid. 

A CSMP program was written (Appendix G) with the above 
equations, which represent the aerodynamic models and the 
control laws of the linear CBTT autopilots. 

Figure 3.7 shows the sideslip angle 8 for the 
linearized CBTT autopilot, having maximum values of -0.46^ 
and 0.26^. 

Figure 3.8 shows the roll angle response which has an 
overshoot of 4 percent and is about the same with the roll 
angle of uncoupled roll autopilot. The above plots verify 


the results, given in the [Ref. 2]. 
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Fig. 3.7  Sideslip angle (8) Ys Time (t) 
Linear CBTT Autopilot (1 RAD Roll 
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Н. RELATIVE STABILITY OF CBTT AUTOPILOTS OF CIRCULAR AND 

ELLIPTICAL AIRFRAMES 

When the coordinating branch gain Kyp was set to zero 
in CBTT autopilot of the circular airframe, the maximum 
sideslip angle increased to 2.8 degrees. Therefore, for 
the circular airframe the coordinating branch is not very 
effective in helping the yaw autopilot to reduce sideslip 
angle. 

TE Кур is increased to unity, as it is for the CBTT 
autopilot of the elliptical airframe, the maximum sideslip 
remains below one degree for 3 seconds but the autopilot is 
unstable. Modifying of the circular airframe will be done 
in the non-linear CBTT, for minimization of the sideslip 
angle. 

Removal of the aerodynamic cross-coupling (i.e С, i 


B 


C. , C ) had no effect on the sideslip angle. This 
Š 


Y R 
showed that the aerodynamic cross-coupling plays an 
indirect role in determining the quality of sideslip 
control by determining the relative stability of the 
coordination branch or the magnitude of the coordination 
gain Кур. The magnitude of sideslip angle is dependent on 


the nulling effects of two parallel paths shown in Figure 


3.9, as given in (Ref. 2]. 
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The contribution of yaw acceleration ny to the maximum 
8 is negligible and therefore neglected. 8 is formed 
mainly by the subtraction of the kinematic paths of a,p and 
the yaw angular rate r. 

The coordination is obtained in the CBTT control law by 
comnanding the yaw autopilot with a yaw angular rate 
command ne (Figure 3.2) of Кур аыр, which forces r to be 
equal to a,.p and therefore nulling as shown in Figure 
3.9. The nulling process will be accomplished more 
efficiently if Кур =1.0. The reason the sideslip is not 
nulled completely is that пу is not zero and r can not 
equal ro over all frequencies. 

The linearized CBTT autopilot of the circular airframe 
was unstable, when the uncoupled channel control laws used 
with a coordinating branch gain Кур of unity. The 
autopilot was stabilized by decreasing to 9.458. The 
actuator of roll control channel has a frequency of 60 
rad/sec. 

The stability of the CBTT autopilot for the elliptical 
airframe is considerably better, due to the setting of Кур 


at unity, which results in a substantial decrease in the 


magnitude of sideslip angle as shown in Figure 3.9. 
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I. CONCLUSIONS 


l. The linear coupled autopilots verify the results of 
the uncoupled autopilots and keep the sideslip angle small 
for coordination motion. (section II.D. and II.E) 

2. Further Conclusions about the relative stability of 
the autopilots, about the gain Кур of the coordination 
branch and about the sensitivity of the autopilots to the 


aerodynamic cross-coupling are given in chapter 7 of [Ref. 


2]. 
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IV. NONLINEAR ANALYSIS OF CBTT AUTOPILOT FOR CIRCULAR AND 


ELLIPTICAL AIRFRAME 


A. GENERAL 

The three dimensional nonlinear aerodynamic models 
used for the analysis are shown in the Figures 4.1 and 4.2. 
For the configuration under consideration, the assumptions, 
which have been found to be consistent are referred to in 
chapter 5 of [Ref. 2]. 

Nonlinear aerodynamics and mass parameters values are 
presented in Appendixes A and B, as taken from tne [Ref. 
2]. The same flight conditions (i.e 60kft altitude and 
M=3.95) as in linear CBTT. 

In this section the control laws for the non-linear 
CBTT autopilots are discussed. The results for commanding 
the CBTT autopilots for both airframes are shown and the 
desired aerodynamic model to enhance CBTT performance is 
determined. The performance of the CBTT autopilot for the 
elliptical airframe is determined, when no lateral 
aerodynamic cross-coupling exists and with ideal airframe 
dynamics. Also the performance of the ?qCBTT autopilot for 
the circular airframe is determined, when no inertial and 
kinematic cross-coupling exist into the pitch autopilot. 


The conclusion of the nonlinear analysis are stated. 
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ig. 4.2 Nonlinear Lateral (Roll/Yaw) channel dynamic 
model 
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Б. CONTROL LAWS FOR ELLIPTICAL AIRFRAME 

The control laws used for the following nonlinear 3-D 
studies were the same as those used for the linear studies 
in section III, which are shown in Figures 4.3 and 4.4 
taken from [Ref. 2], except for the gain as shown in the 
bold line of the coordination branch in Figure 4.4. The 
new gain a is held constant at one degree magnitude for 
angles of attack less than one degree positive and greater 
than negative 5 degrees. For angles-of-attack greater than 
one degree positive, the gain a is equal to the angle of 
attack. This maintains coordination for very small 
angles-of-attack. 

Gravity effects were not included in the linear 
studies, because it was assumed to nave a negligible 
influence on autopilot stability and response perturbations 
about a missile trim conditon. 

Gravity effects were included in the nonlinear 
studies, where the missile body-fixed yaw axis will be 
subjected to the full force of gravity and may therefore 
have a significant influence on sideslip. 

1. Pitch Control Law 

The pitch control law for the elliptical airframe 
is the same that is used in the linear studies in section 


III.F., except in the command applied to the channel. 
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In inertial rectangular coordinates, the 
acceleration commands, are: 


а. пе (acceleration command in inertial Zy 


direction as in figure 2.4) =n, - cos? mye aa 
С 
where: n, : guidance command (gees) 
с 


-cos? : anti-gravity bias command (gees) 


9: Elevation Euler angle =  (qcos?-rsin?)dt 


bon acceleration command in inertial 2 


B 
C 


direction, which is equal to zero, because it is desired 
coordinated motion. In the polar coordinates, the 
acceleration commands are: 

- | (m, - cose) + ny? ] É 


– ( Ст, - 56) 1% 
- (Na. - со50) 


а. Vic 


(1057152) 


р Фес tam (©) 
ic (IV.B.1-3) 


hence 2 cong andenoso.i2..- 

Since the pitch control law of the elliptical 
airframe does not have an integration in the acceleration 
error path (i.e in the acceleration compensator as in the 
case of circular airframe) and requires a gain in series 
with the acceleration command shown in Figures 4.3 it was 
necessary to modify the antigravity command as follows to 


assure an anti-gravity bias of just one gee. 


T21 





n (elliptical airframe inertial acceleration command) = 


= n, - 0.913 cose (LIT B 1.4) 
С 


The autopilot was tested, using a commanded 2 gee 
(0°, 180°), that denotes an inertial command of 2 gees, 
which is first applied in the ^ ог upward direction at T = 
2 seconds for 3 seconds and at the fifth second is applied 
in the 180° or downward direction. Since at first both the 
missile roll angle and roll angle command are at zero 
degrees, there is no roll motion and the missile turns 
upward as a skid-to-turn controlled missile. At fifth 
second the missile is commanded to roll through 188 degrees 
while moving in a coordinated manner in yaw and roll to 
minimize sideslip angle and prevent or minimize negative 
angle-of-attack. 

The other equations, which represent the nonlinear 
pitch control law are the same as for linear studies 
section III.C. and the state variables are shown in Figure 
з. 


The equations аге: 


X 2 -150X + 150 N25 (IV.B.1-1) 
Ү = -6ү- 0.5305 1, + 0.13 X (ТОБЕ 1=2) 
р. = 0.383 (У+9 151:3)-3.0+(ү- 915.5) (1У.в.1-3) 
Sp = - 1984 др + 10395.3 др, ЕТЕ) 
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where: 


achieved body-fixed acceleration in Z axis 


n : inertial acceleration command. 


2. Lateral Control Laws 

The lateral control laws are the same as those 
used for the linear studies, except the modification in the 
coordination branch. 

The command ee for the roll control law is zero 
for the first 5 seconds and 1809 affterwards, as derived in 
wv.B.1l. 

The equations which represent the lateral control 


laws are the same those of section III.F. 


e 


Yo = -4 Y 99: Муд (IV.B.2-1) 
Е (ae- py “ _ аеР 
ду, = 0.608 ( Y, E E ө: ted 53 UNS. 
- Y M а : 
= 0.603 (а. = 33 e o => - 
(IV.B.2-2) 
where: 
de = о | a> ° 
B cu 
and 
Qe = 4 a < | 
Oe = О 
Х, = - а, - њеф+ 116 Фе AS 
dy = -198.4 SY +10395.3 дус MN 
Yı = 5% +0.05035Cx,-Dis33) +0455 (X4-Plsr.3) ER 
Nc -6X, * 0.018 (P[51.3) mS 


др = ~ 15 дес + 10425 (W -X2)462.55 (Yi-X2) (ТУ.В.2-7) 





бр = –188.4.др + 10%95,3 дес VEI 2-01 


The state variables are shown in Figure 4.4. For 
purposes of analysis, a CSMP program was written (Appendix 


H). 


om AERODYNAMIC MODELS FOR ELLIPTICAL AIRFRAME 
The aerodynamic models of Pitch and Roll-Yaw channels 
for the elliptical airframe are shown in Figures 4.1,4.2. 
The airframe dynamics under consideration, have a level of 
complexity sufficient to determine the critical areas of 
concern regarding the stability and response of CBTT 
E»ntrol. 
1. Pitch Channel 
The nonlinear aerodynamic model of the Pitch 
channel consisted of the uncoupled model, coupled with the 
Roll-Yaw channels via the inertial cross-coupling of rp 
into q and the kinematic cross-coupling of -8.P into a. 
Also the antigravity command cos% cos9 is inserted into n 
This model is shown in Figure 4.1. 
The above terms consist of the nonlinear terms of 
the pitch channel. 
The equations which represent the pitch 


aerodynamic model are: 
q P.r|st3+ 51.3 98d Cu [Toy poss 
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À =Q +K. g — P.6/57,3 (IV.C. 1-3) 


where: 
1845 
K = — = 0.48 
V 
а = 1659 (1ЬЕ/ ЕЕ.) 
5 = т (Et?) 
d = 2 (ft) 
W = 2475 (lbf) 
2 
i = 790 lug- 
уу (slug-ft 


Functions Cu Gr 3g) and Canlar fp are nonlinear 
functions shown in Appendix B that vary with a and ӛр. 
Substituting the values into the above equations, 


they yield: 


N, = —1.3333 - (y (IV.C. l-4) 
д = 0.19452: Р-ү + #51. 9541 Си (IV.C.1-5) 
aa = Q - 0. 659438 CN - 0.014452 P- € (IV.C.1-6) 


2. Lateral Channels 
The nonlinear aerodynamic model consisted of the 
uncoupled roll and yaw channels, which coupled via C. 
B 


en and Cn ». Also the lateral channel coupled with 
6 6 


Y R 
the lateral channel coupled with the pitch channel via the 
inertial cross-coupling of -qp into r and the kinematic 
cross-coupling of ap into 8. The above term is the 


nonlinear part the lateral pitch channel. The nonlinear 


lateral dynamic model is shown in Figure 4.2. The 
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equations that represent the aforementioned model are as 


follows: 
We = 4 5 (6: + Cygy -dy) | W (gees) (IV.C.2-1) 
8 = Km + aplsta -r (deglsec) (IV.C.2-2) 
+ = -9-Р1[5%3 +545 454 (сш. +, + $ + Cog, Se) [Tez 
(о 5) 
Р = 55954 (Сасе + G би + Су -5е) [1 ау.с.2-4) 
Ф = P/srt3 (xadlsec) (IV.C.2-4a) 
Ф, = Р Cdeg (sec ) (Fu re oss) 
where: 
K = 0.48 
а = 1650 (1bf/ft,) 
S =n (£t?) 
а = 2 (ft) 
w = 2475 (1bf) 
I, = 110 (slug-ft* 
БЕ - 953 (slug-£t7) 


as referred to in Appendix A. 


Functions C (a), m: С, (а), С, (а), с, 


Y 
В бұ бұ 


Š y 


С, (а), с, (а) апа С (а) are nonlinear functions which 
n 
8 5n 5n 


vary with a angle-of-attack and referred in Appendix B. 


Substituting the values into the referred before 


equations one obtains: 
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мүр = 4.3333 ( €. 8 + Суу. у) tia mD 


B = 0.63998 (04.8 36, бу) “0.0452 а-Р-Ү/ (1У.с.2-6) 
( =-0.0845 4: Р+69642 (О, 6 +Су dy +Су, бе) (1у.С.2-7) 


р = 5400.397 (Gg. E Cg y € Cop, 08) АЯ) 
Ф = P Ие 2209) 


D. ANALYSIS OF NONLINEAR CBTT AUTOPILOTS FOR ELLIPTICAL 

AIRFRAME 

For purposes of analysis, a CSMP program was written 
(Appendix H) with the aforementioned equations of the 
aerodynamics models and control laws of the nonlinear CBTT 
autopilots. 

To decrease computer time, appropriate initial 
conditions were inserted. The initial conditions are: 

a : angles of attack = 2.41 degrees 
6р : pitch tail incidence = 0.658 degrees 
8 : pitch Euler angle = 0.636 rad 


output of pitch acceleration feedback laq = lgee 


>< 


Y : output of pitch acceleration compensator = -0.0105 


8 : pitch actuator command = 0.0148 rad 


Pc 


The achieved maneuver plane acceleration п. is 


calculated from the body-fixed accelerations De (Figure 
В 
4.6) and (Figure 4.7) and the roll angle y (Figure 4.13) 


as follows: 


r27 





NZ (GEES) 


Fig. 4.5 Achieved Inertial Acceleration (n 
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CBTT of Elliptical airframe: 2 Geés (0°, 18@ 





) vs Time 


t) 
5, 
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Fig. 4.6 Achieved Body-fixed Acceleration (n, ) vs 


Time (t) CBTT of Elliptical airframe; 2 Gees 
(07,1807) 
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NYBIGEES) 





Fig. 4.7 Achieved Body-fixed Acceleration (ny ) vs 
B 


Time (t) CBTT of Elliptical airframe; 2 
Gees (Ø , 188°) 
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М> = N, cosb + Мү. • sim $ (IV.D-1) 


Figure 4.5 through 4.16 show the responses of noe 


ny , As В, а, р, К, Š, бр, б, апа бр. Also these figures 
% table II show the time constants and the overshoots of 
the referred before variables. 

Figure 4.5 shows that the achieved inertial 


acceleration n, doesn't satisfy the requirement of the 


2 
overshoot (i.e below 18 percent) and it is not acceptable. 
The results of this analysis verify the result of the [Ref. 


2]. 


E. CBTT PERFORMANCE WITH IDEAL AIRFRAME DYNAMICS FOR 

ELLIPTICAL AIRFRAME 

The purpose of the following simplifications to the 
airframe dynamics model is to isolate the critical 
cross-coupling paths, which have caused, the transients in 
the maneuver plane acceleration responses. 

The dynamic model without the coupling paths will be 
referred to as ideal dynamics. Although ideal dynamics are 
not physically attainable, it is a useful goal for both 
autopilot and airframe designers. 

The same guidance commands are applied to the CBTT 
autopilot of the elliptical airframe as in par. ІУ.В.1 
(i.e, 2 gees (9°, 180°) but with lateral aerodynamic 


cross-coupling removed. That means, the aerodynamic 


TL 








Fig. 4.8 Angle of Attack (a) vs Time (t CBTT_ of 
Elliptical airframe; 2 Gees (Ø`, 189°) 
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Fig. 4.9 Sideslip angle (g) vs Time (t) 
re of Elliptical airframes; 2 Gees 
(0 189° ) 
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Fig. 4.10 Body pitch Angular rate (а) vs Time (t) 
CBTT of Elliptical airframe; 2 Gees 
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Fig. 4.11 Roll angular rate (p) vs Time (t) 


CBTT of Elliptical airframe; 2 Gees 
(09, 1809) 
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Fig. 4.12 Yaw angular rate (r) vs Time (t) 
| CBTT of Elliptical airframe; 2 Gees (09. 180°) 
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Fig. 4.13 Roll angle (4) vs Time (t) 
CBTT of Elliptical airframe; 2 Gees (9%, 1897) 


Jo 

















Pitch tail incidence ($)) vs Time (t) 


CBTT of Elliptical airframe; 2 Gees (0°, 180°) 
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| Fig. 4.15 Yaw tail incidence (бу) vs Time (t) 
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CBTT of Flliptical airfframe; 2 Gees (Geer 1807) 
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Fig. 4.16 Roll tail incidence ($4) vs Time (t) 


CBTT of Elliptical airframe; 2 Gees (42, 189^) 
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cross-coupling was removed. The major influence of the 
aerodynamic cross-coupling in the elliptical airframe 
response has been to decrease the overshoot in the achieved 
maneuver plane acceleration resulting from the second 
guidance command. 

In addition to aerodynamic cross-coupling the 
kinematic cross-coupling of -8p into a and inertial 
cross-coupling of pr into q were removed. Therefore, the 
only cross-coupling which exist in the airframe dynamic 
model are the kinematic coupling of ap into 8 and the 
inertial coupling of -qp into r. There is also the 
autopilot cross-coupling of the coordinating command from 
the Roll to Yaw channel. 

The equations, which represent the modified 
aerodynamic model by removing the referred cross-coupling 
are same with those aS in paragraphs IV.B.1, IV.B.2, ІУ.С.1 
and IV.C.2, with the only difference in the following 
equations: - 

1. Pitch Aerodynamic Model 

Q = 51.3 Q Sd Cu [Tyy (Чед [52<7) 
= 151.9541 Си (IV.E.1-1) 


and e 


a =- 9+ Kg 
= q- ( K á 8 lw) 


= 9- 0.65998 Су (degísec ) (IV.E.1-2) 
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2. Lateral Aerodynamic Model 
ес- 4.Р 4 57.3 4,54 ( az. 8 + (ту ду) 


61.3 I22 
= -004452 q- P 4 696.43 ((ng- 8 * nz, Y) (deglsec) 
(IV. E 1=3) 
and s 
· 5:38 "esa 
р = Туу Се, др 
= 5400.3377 Cep -ðr (deglsec”) 
Ё ПӘЛЕ 2-2) 


For purposes of analysis, a CSMP program (Appendix 
H) was written, using these equations. Figures 4.17 
through 4.24 show the time responses of the n 


n а, B, 


2! Z5 
E Y, бр, бу, бр. The above figures and table II show the 
time constants and the overshoots of the variables. 

Figure 4.17 and Table II show that the achieved 
inertial maneuver acceleration satisfy the requirements for 
the time constants and overshoot, when the cross-couplings 
which were discussed before have been removed. 

Figure 4.25 (taken from [Ref. 2] shows the 
critical feedback paths which couple the pitch and yaw 
channels via missile roll rate. 

The aforementioned figures and plots verify the 


results of the [Ref. 2]. 


P. CONTROL LAWS FOR CIRCULAR AIRFRAME 
The control laws used for the circular airframe are 


the same as those used in the linear studies, except for 
142 
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Fig. 4.17 Achieved Inertial Acceleration (nz) vs Time (t) 


CBTT of Elliptical airframe; Ideal 
Airframe dynamics; 2 Gees (99,189?) 
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Fig. 4.18 Angle of Attack (a) vs Time (t) 
CBTT of Elliptical airframe; Ideal. 
airframe dynamics; 2 Gees (Ø`, 189°) 
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Fig. 4.19 Sideslip angle (3) vs Time (t) 
CBTT of Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (0°, 18@°) 
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Body pitch angular rate (q) vs Time (t) 
eB or Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (8 , 180°) 
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. 4.21 Body Yaw angular rate (r) vs Time (t) 


CBTT Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (8%, 180°) 
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Fig. 4.22 Pitch tail incidence ($ ) vs Time (t) 
CBTT of Elliptical airffame Ideal 
airframe dynamics; 2 Gees (0°, 189°) 
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Fig. 4.23 Yaw tail incidence (8_) vs Time (t) 
CETT of Elliptical airframe; Ideal 
airframe dynamics; 2 Gees (02, 1809) 
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Fig. 4.24 Roll tail incidence (6,) vs Time (t) 


CBTT of Elliptical airframe; Ideal 
i ' о 
airframe dynamics; 2 Gees (0 , 180 ) 





Pig. 
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4.25 





Critical Feedback loop for CBTT 
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TABLE II. CBTT PERFORMANCE OF ELLIPTICAL AIRFRAME 


VARIABLE VALUE TIME (sec) CROSS- T, (sec) (sec) 
COUPLING  OVershoot$ Overshoot’ 

nomin -2.934(gees) 5.061 A11 0.48 and 0.41 and 

К тах 1.351 (qees) 5.661 A11 0% 11.25% 

nomin -2.934(gees) 5.062 ideal 9.48 and 9.37 and 

пах 9.937 (дее5) 5.893 ideal 0% 0% 

As 1.35(deg) 6.10 All 

EN 6.75 (deq) 3.18 All 

BE in 2.22(deq) 5.295 ideal 

EE ax 6.75 (дед) 3.18 ideal 

in -8.9(deg) 5.49 All 

Bax 1.49 (deg) 6.05 All 

Bein 1.53(deg) 5.45 ideal 

Е. 1.87 (deg) 5.91 ideal 

dein -18.6(deg/sec) 5.639 All 

Bux 9.5(deg/sec) 25337 All 

gin -7.8(deg/sec) 5.328 ideal 

ШЕ, 9.5(deg/sec) 20337 ideal 

in -14.115(deg/sec) 5.600 All 

— 282.094 (deg/sec) 5.387 All 

in ЕШ 25(464/ ес) 6.600 ideal 

m 281.928(deg/sec) 5.387 ideal 

Din -2.32(deg/sec) 5.99 A11 

Ex 33. 1 (деа/5ес) 3.032 All 

TN -0.178(deg/sec) 7.036 ideal 

Dx 37.5 (deg/sec) 5.29 ideal 

Note: 


All the above referred values are for command 300, 180°) 


тү: 62 percent time constant of achieved maneuver plane 
acceleration due to the first guidance command. 

тә: 63 percent time constant of achieved maneuver plane 
acceleration due to the second guidance command. 
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the gain (a) of the coordination branch of lateral controi 
laws, as mentioned in the study for the elliptical 
airframe. 

The same 2 gees (9°, 188°) guidance commands and 
flight condition were used for direct comparison with the 
performance results of the elliptical airframe. 

However, to decrease simulation run time, the first 
guidance command and anti-gravity bias were applied at zero 
time with no missile or autopilot initial conditions. 

1. Pitch Control Law 

The pitch control law for the circular airframe is 
the same as that use in the uncoupled autopilot shown in 
Figure 4.26. The command used in the autopilot is 
different than the one used in uncoupled channel and it is 
as follows: 

Ус = Иа. = (950 СЕРЕГЕ И) 
as referred to in section ІУ.В.1 

The equations which represent the nonlinear pitch 
control law are the same as those in section III.C. of 
the linear pitch control law, using Ka, -0.0387 


The equations are: 


x z-150 X -- 150 2g (IV.F.1-1) 
Y= Yy SS (ме-Х) (IV.F.1-2) 
Yu 2 -SYy + 6.572-(Nc-X) (IV.F.1-3) 


др,= - 0445. бр, = 0. 3u058( Y, - à [51.3) -2.2308-(У3- 9193) 
(TV. P.1—-4) 
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др ш -188.Ң Әр + 10795.3 „фр, (IV.F.1-5) 


The state variables of the above equation are 


shown in Figure 4.26. 


2. Lateral Control Laws 


The lateral control laws for the circular airframe 


are the same as those used for the linear studies, except 


the modification in a 


airframe. 


(or a) which was referred to before. 


The same command is used, as in the elliptical 


The equations which represent the lateral 


control laws are the same those of section III.C.1 and 


ШІ.С. 2. 


Y, 2 -5% + 16 Wya 


and they are: 


(IV.F.2-]) 


бу, 0.435 sl TASE . (de P+ Xe P )] +35 (17+ 93-28082) 


where: 


and 


573 (5737 


= 0.485 [Ү; + E - 0.000433 (dp P+ de: 2] +435 (+5 -о. 00139 «e-P ) 


(IV.F.2-2) 

Се = < 6 

E^ | к>! 

Zza 

del $ к<! 

Qe =O 
ду = - 188.4 dy + 10795.3 dy, (IV.F. 2-3) 
XY, 2 -8X -116 16 d (IV. F.2-4) 
Y, = -5 Y, + 0.05033 (x, = P/s7.3) 10.155 (X -P/57.3) 


(IV. Е. 2-5) 


155 





X2 2 -6X2 t O.012 ( P/51.3) (IV.F.2-6) 
SR = 15 de -- 0.25 (Y,- X2) 4 45 4 -X2) (IV.F.2-7) 
др = -19$46р. +40495.3 -дре (IV.F.2-8) 


The above state variables are shown in Figure 


G. AERODYNAMIC MODELS FOR CIRCULAR AIRFRAME 

The aerodynamic models of Pitch and Roll-Yaw channels 
for the circular airframe are the same as those of the 
elliptical airframe and are shown in Figures 4.1, 4.2, the 
only differences which exist between them, are the 
aerodynamic mass parameters and nonlinear derivatives 
changes, as referred to in the (Appendices A and B). 

Hence, the equations which represent the dynamic 
models are: 


1. Pitch Channel 


Vag =-QS Cul w (gees) (ІУ.С.1-1) 
a P. Y 533 3 59 Си де se 2 (IV.C.1-2) 
a + se (degisec 

P.6 


a=Q + И. Изв — (deg |sec) (IV.C.1-3) 


53.3 
Substituting the | for aerodynamic and mass 
parameter for circular airframe, the above equations 
become: 
Yap = = o 05252 - C (IV.G.1-1) 
Q Z-0.047452 PY+ 436.3604 Cu (IV.G.1-2) 
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& = q - 0.9854 Cy- 0.017452 P- (IV.G.1-3) 


2. Lateral Channels 


W, = 98 (Су. + Or, -öy) Iw (gees) (пе =) 
4 - K-Nyg ROS PÍSES - Y (deglsec) 
v + ( Gag: Саду: бу + Eng, $2) 
ES 2 5, 
Б = 513 35d | | Ce, - 6 4 Gp, dv + Ces, -öR) 
Ixx 

м. с. 2-4) 
Фа ES P (deg [sec ) (IV.C.2-5) 

Substituting the values in the above equations, 

they yield: 

Wa = 2.05292 ( 4.6 * Os, .5у) (IV.G.2-1) 


$ = 0.9954 (CyB + Gy dy) +0075 Pa-F (1У.6.2-2) 
ү =-0.013452 Ч-Р +353.385+ (Сл: 6 + Спр-бу +Саб; -92.) 
УС. 2-3) 


H. ANALYSIS OF NONLINEAR CBTT AUTOPILOTS FOR CIRCULAR 

AIRFRAME 

For purposes of analysis, a CSMP program was written 
(Appendix H) with the equations representing the 
aerodynamic models and control laws of the nonlinear CBTT 
autopilots for referred airframe. 

At first, in order to understand the effect of faster 
response, in the reduction of the effect of kinematic and 


inertial cross-coupling, the pitch acceleration error gain 
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мана 





Ky Figure 4.26 was reduced to K, = 0.9274 for slower 
1 
response. 


The equation, which represents the acceleration 


compensator after the above reduction, is: 


ya - 10:0274 ($10.143 — 1) 
Ë S 
po) 


Utilizing method of state-space representation with 


e (4-9 ГЕЛ» Не) 


forcing term derivatives, the equation IV.H-l yields: 


Ya = Yu - 0. 95 804 - (Nc-X) (IV.H-2) 
Yy = -5-Y +4.653 MX) НЭ) 


Using the above equations, instead of the equation 
IV.F.1-2 and IV.F.1-3, a new CSMP program (Appendix I) was 
written. 

Figure 4.28 shows that the achieved maneuver plane 
acceleration has a good response to the first command which 
is applied at zero time. The missile with the circular 
airframe moves upward like a skid-to-turn missile, as did 
the elliptical airframe missile, because the motion is in 
the deSired maneuver direction and therefore the roll 
channel is not commanded. 

The time constant of the achieved maneuver response ті 
is equal to 9.49 seconds. The achieved maneuver plane 
response due to the second guidance command applied at 3 
seconds, shown in Figure 4.28 is reacting differently to 


the kinematic and inertial cross-coupling than the 


T59 
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-1.0 





0.0 0.5 1.0 15 20 25 30 3.5 
TIME(SEC) 


Fig. 4.28 Achieved Inertial Acceleration (п) vs Time (t) 


4.9 4.5 5.3 3.3 6.0 


CBTT of Circular airframe; 2 Gees ІШЕ! 180°), 


ña 
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elliptical airframe in Figure 4.5. Rather than overshoots 
and undershoots, a slowing transient starts at 3.5 seconds. 
Figure 4.29 shows that the overshoot in the body-fixed 
pitch acceleration, due to the kinematic and inertial 
cross-coupling during the second command, is 698.8 percent 
which is substantially more than it was for the elliptical 
airframe and occurs much sooner at 3.66 second. 
Figure 4.29 through 4.39 show the time responses of 
the ny , a, 8, дур, г, %, бр, Š and Sn. Above figures 
and En. III show the characteristic of the responses of 
the variables for the circular airframe, when the K, = 
-0.0274. 


The referred results of this study verify the results 


of [Ref. 2]. 


Í. EFFECT OF INCREASING PITCH CHANNEL SPEED OF RESPONSE 
In order to reduce the effect of kinematic and 

inertial cross-coupling during the second guidance command, 

the response of the Pitch channel of the CBTT autopilot for 

the circular airframe was made faster as shown in section 

II. This was accomplished by increasing 

the acceleration error gain K, - -0.0387. 


2 


Using the above gain K the equations represent the 


А 
2 
acceleration compensator are the IV.F.1-2 and IV.F.1-3. 


The effect of the change in K, on achieved body fixed 


acceleration (Figure 4.40) results in the achieved maneuver 
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0.0 05 1.0 15 20 25 30 35 40 45 50 5.5 6.0 
TIME(SEC) 


Fig. 4.29 Achieved Body-fixed Acceleration (n, ) vs 
B 


Time (t) CBTT of Circular airframe; 2 Gees (9°, 
1809), К. = -0.0274 
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' €.0 OS 10 15 20 25 3.5 40 45 50 35 6.0 
TIMEISECH 


Fig. 4.30 Achieved Body-fixed Acceleration (ny ) vs 
B 


Time (t); CBTT of Circular airframe; 2 Gees (0°, 
180%), K, = -0.0274 
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Fig. 4.31 





15 20 25 30 35 40 45 $0 5.5 
TINELSEC) 


Angle of attack (a) vs Time (t) 
CBTT of Circular airframe; 


2 Goes (0, 180°), КА = -0.0274 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 5.0 


ТІМЕІ5ЕС! 


Fig. 4.32 Sideslip angle (g) vs Time (t) 
CBTT of Circular airframe; 
2 Gees (0%, 180%), K, = -0.0274 
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2.5 3.5 40 4.5 8.0 
TIME (SEC) 


-10.0 


Fig. 4.33 Pitch Body angular rate (q) vs Time (t) 
CBTT of Circular airframe; 2 Gees 
(09, 180 ), К, = -9.0274 
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Fig. 4.34 Yaw Body angular rate (r) vs Time (t) 
CBTT of Circular airframe; 


28cces (0, 180), КА = -0.0274 
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Big. 4.35 Roll Body angular rate (P) vs Time (t) 


CBTT of Circular airframes; 
2 Gees (@°, 180%), K, = -9.0274 
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Fig. 4.36 Roll angle (4) vs Time (t) 


ОВТТ ОЕ Circular airframe; 
2 Gees (9° ; 189° ); КА = -0.0274 
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Fig. 4.37 Pitch i incidence (6 >) vs Time (t) 


CBTT of Circular airframes; 
2 Gees (Ø, 180), K, = -0.0274 
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Fig. 4.38 Yaw tail incidence (6) vs Time (t) 
CBTT Circular airframe; 


2 беев (09, 1809), КА --0.0274 
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Fig. 4.39 Roll tail incidence (p) vs Time (t) 
CBTT Circular airframe; 
2 беев (09, 1809), Кл --й0.0274 
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TABLE III. 


VARIABLE VALUE 


Nzmax 
Amin 
п2 пах 
Nzmin 
П7тах 
min 


a 
max 


a 
min 


a 
max 


a 
min 


a 
max 


a 
min 


B 
max 


Bmin 


max 


max 


max 
min 
max 
min 
max 
min 
max 
min 
max 


min 


1.092(gees) 
-3.844 (gees) 
1.189 (gees) 
-3.125(gees) 
1.882 (gees) 
=3. 125 (gees) 
11.15(deg) 

ё (дед) 
11.489 (deg) 
8.853(deg) 
11.488 (deg) 
3( deg) 
4.25(deg) 


-8.399 (deg) 


3.633 (deg) 
-8.326 (deg) 
5.114 (дед) 

08.414 (дед) 
21.39 (deg/sec) 
-8.84(deg/sec) 
29.562 (deg/sec) 
-14.323(deg/sec) 
29.562 (deg/sec) 
-17.575(deg/sec) 
278.5 (deg/sec) 
-14.166(deg/sec) 
279.487 (deg/sec) 
-14.871(deg/sec) 
277.574 (deg/sec) 
-13.798(deg/sec) 
22.66 (deg/sec) 
-1.456(deg/sec) 
21.486 (deg/sec) 
-1.322(deg/sec) 
23.159 (deg/sec) 


-1.364(deg/sec) 


CBTT PERFORMANCE OF CIRCULAR AIRFRAME 


TIME(sec) 


4.548 
8.825 
5.848 
8.825 


1.432 


8.862 
3.656 
8.862 
8 

3.461 
4.286 
3.458 
4.298 
3.688 
4.628 
6.285 
3.228 
0,261 
3.225 
8.261 
3.251 


3.394 


CONDITION 
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v, (sec); 
overshoot’ 


8.46 and 
% 
0.34 and 

5% 
0.34 апа 


5% 


(Sec); 
overshoot? 


8.52 and 
% 
8.62 and 

0.958 
8.53 and 


0% 









0.0 0.5 









1.0 1.5 20 25 30 35 40 45 50 55 6.0 
_TIMELSED) 


паса [ner ial Acceleration (n,) vs — 
Time (t) CBTT of Circular airframe with faster 


responding Pitch channel; K, = —0.0387, 2 Gees 
(047, 180 ) 2 
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plane acceleration response (Figure 4.41). Figure 4.41 


shows that the achieved maneuver plane acceleration 
response during the first guidance command has improved. 


Table III shows that increasing the K, resulted in a 
decreasing in maximum sideslip angle 3.66 and a slight 


decrease in yaw angular rates (maximum 21.4 degrees). 


EFFECT OF INERTIAL AND KINEMATIC CROSS-COUPLING IN 


J. 


PITCH CHANNEL 
intosg , 


The cross-couplings, (~8P) into a and (rP) 
were removed in order to assess their effect on 


performance. Hence, the equations for a and q in the pitch 


aerodynamic model becomes: 
EU 


% = q + x. Mg (deg [sec) 


Replacing the values to the parameters of the above 


(IV.J-2) 


equations, they yield: 
а = 9-0. 9354 (м (1У.2-3) 
(ІУ.2-4) 


Q - 138.860 Cu 


nonlinear functions which vary with a and 


where Сы, Cn 
The lateral aerodynamic 


Š referred to in (Appendix B). 
For purpose of analysis, a 


|Р 

cross-coupling was retained. 
CSMP program (Appendix I) was written. 
Comparing Figures 4.42 and 4.43 with figures 4.40 and 


4.41 it is concluded that the undersirable transients аге 


gone. 
175 
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. 4.41 Achieved DES Acceleration (n, ) vs Time 
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(t) CBTT of Circular airframe; Faster responding 
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Pitch channel; 2 Gees (0', 180 ) К, = -0.0387 
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“TIME (SEC) 


4.42 Achieved Inertial Acceleration (nz) vs Time (t) 


CBTT of Circular airframe; Faster Pitch 


channel; inertial and kinematic cross- ‚coupling 
into Pitch channel removed; 2 Gees (09 i 1809 IR 


ae -0.0387 
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. 4.43 Е Body- "n Acceleration (n, ) vs Time 
B 


(t); CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross- coupling 
into Pitch channel removed; 2 Gees (а ој y 
КА = - 0. 0387 
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Table III shows that the achieved maneuver plane 
acceleration time constant for the second guidance command 
has decreased. 

Figures 4.44 through 4.50 show the time respanses for 
ШО, а „т, бр, бу, апа бр when the Кл is equal KAS -0.0387 
and the kinematic and inertial cross-coupling of the Pitch 


channel have removed. 


К. CONCLUSIONS 

1. The maneuver plane responses have transients, 
which may have to be reduced, are caused by inertial and 
kinematic coupling between pitch and yaw dynamics. 
Transients cause excessive overshoots and undershoots in 
achieved maneuver plane acceleration for elliptical 
airframe (Figure 4.6) and excessive slowdown in the speed 
of response for circular airframe (Figure 4.28, 4.40). 

2. The result of the nonlinear 3-D performance study 
verify the linear study of chapter III. 

3. The results of this study verify the results of 
mer. 2]. 

Further analysis and conclusion about a, q, r, P 


and 8 are discussed in [Ref. 2]. 
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. 4.44 Achieved Body-fixed Acceleration (пу ) vs Time 
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(t) CBTT of Circular airframe; Faster pitch 
channel; inertial and kinematic cross- re 
into pitch channel removed; 2 Gees (60 130°); 


КА, -ø . 0387 
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Angle of Attack (a) vs Time (t) 

CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees (0°, 180°), 
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Fig. 4.46 Sideslip angle (g) vs Time (t) 


CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross- 
coupling into Pitch channel removed; 


2 Gees (6°,188°%), K, = -0.0387 
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Pitch Body angular rate (q) vs Time (t) 
CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees м 189°), 
K, = —0.0387 
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Yaw Body angular rate (r) vs Time (t) 

CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees (0', 1899), 
К, = -0.0387 
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Pitch tail incidence ( 85) vs Time (t) 


CBTT of Circular airframe; Faster Pitch 
channel; inertial and kinematic cross- coupling 
into Pitch channel removed; 2 Gees (9° 7 189° E 
K, = -0.0387 
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Yaw tail incidence ( &y) vs Time (t) 


CBTT of Circular airframe; Faster Pitch 


channel; inertial and kinematic cross-coupling 
into Pitch channel removed; 2 Gees (9°, 180^), 
К = -0.0387 
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V. MINIMIZATION OF THE KINEMATIC AND INERTIAL 


CROSS-COUPLING EFFECTS 


A. GENERAL 

The nonlinear CBTT autopilot for both airframes, 
consisted of the uncoupled autopilots discussed in section 
II, which coupled via inertial, kinematic and aerodynamic 
couplings. 

The effects of the cross-coupling were discussed for 
both airframes in section IV. 

In this section the kinematic, inertial and 
aerodynamic coupling effects will be sumarized. The 
kinematic and inertial effects in the pitch channel will be 
analyzed. 

The minimization of the above coupling and their 
effects will be studied, using feedbacks of angle of attack 


and rate of angle of attacks. 


B, CROSS-COUPLING OF CBTT AUTOPILOTS 

As discussed in previous sections, the four 
cross-couplings, between the autopilot channels for both 
airframes are: 
1. Aerodynamic Cross-Coupling 

The aerodynamic cross-coupling is the coupling 


between the Roll and Yaw nonlinear channels, via the 
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Су г С, апа en (Figure 4.2). The major effect of the 


Š Š 
: Y R 


above coupling in the elliptical airframe response is to 
decrease the overshoot in the achieved maneuver plane 
acceleration resulting from the second guidance command, 
(i.e 2 gee (1899), the other effects on the missile 
variables were small but in a direction which improves 
performance (i.e decreased missile body angular rates, less 
control surface motion, less sideslip variations). In the 


circular airframe, the effects of the C, , C, are 
8 Š 
R 


negligible, the critical coupling coefficient is the C, 
Š 
which effects the stability in the coordination branch, 


which can be improved by reducing the effect o£ C, on 


Š 
Y 


missile roll angular acceleration. (equation IV.C.2-4) 
2. Cross-Coupling Between Roll-Yaw Autopilots for 

Coordination Motion 

This cross-coupling branch, shown in Figure 4.4, 
has been added between the yaw and roll control laws, in 
order to provide coordinated motion between these two 
channels. For purposes of stability the coordinating 
anch gain Кур was selected for the elliptical airframe 
Kyp = 1 and for the circular airframe Кур = 0.458. The 


elliptical airframe is more stable than the circular with 


the above selected Кур, 
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3. Kinematic and Inertial Cross-Coupling 


The kinematic and inertial cross-coupling exist 
between the pitch and lateral nonlinear aerodynamic models 
of the CBTT autopilots. There exist the following 
cross-coupling paths between pitch and lateral channel. 

a. Kinematic cross-coupling path of (-8p) into a 
in the pitch aerodynamic model (Figure 4.3). 

b. Kinematic cross-coupling path of (ap) into 8 
in the lateral aerodynamic model (Figure 4.4). 

с. Inertial cross-coupling path of (pr) into q in 
the lateral dynamic model (Figure 4.4). 

а. Inertial cross-coupling path of(-ap) into r in 
the lateral dynamic model (Figure 4.4). 

The efffects of the inertial and kinematic 
coupling are the transients in the achieved maneuver plane 
acceleration response, and have the form of overshoots and 
undershoots discussed in section III, IV. Transients may 
also result in a slower maneuver plane acceleration 
response and are less pronounced at higher acceleration 
levels. 

Coupling transients may be reduced by increasing 
pitch stability, by decreasing maximum roll rate (section 
IV), by improving the autopilot coordination technique to 
minimize sideslip rate 8 (section IV), by increasing the 


effects of stabilizing lateral aerodynamic coupling 
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(negative C, ), by increasing pitch channel speed of 
B 
response (section IV) and by increasing the pitch stability 


Са with feedbacks of both angle-of-attack and rate 
a 


of angle-of-attack for the pitch autopilot. 
C. KINEMATIC AND INERTIAL CROSS-COUPLING IN PITCH 
AUTOPILOT 
The kinematic and inertial cross-couplings in the 
pitch autopilot, are via the pitch angular rate q and the 
rate of angle-of-attack (a) (Figure 4.1)). 


The equations for q and a (figure 27) are: 





-_ 51.349 Sd es a УР (V.C.1) 
PF т о (саф P+ (ma; үк SE 
ал 0- (P.8/53.3) (ЕЕ 20 
where: 
rp/57.3 = Inertia Cross-Coupling 
-8P/57.3 = Kinematic Cross-Coupling 


Removing Са from the parenthesis in (V.C.1) and 
a 
combining with (V.C.2) to eliminate q, the (V.C.1) results 


in: 
ee 


_ 543954 cP _ (Pé Y 
E nyn с ы (ma) + s (S) 


5154 54 Спа (2 ао) v.P “сз 
T ma 51.3 57.3 


Ve) 
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Let: 





и. = 54.5 453 | oe 
Iyy 
and ( 
Kr = ge 
Сис. 


Taking the Laplace transform of (V.C.3) and solving 


for a, the equation becomes: 


P.6 
-1 Pls7.3) - S: (3533 
сөл En... (ерја 3) - yy gs UU SG) 
BÉ. 1 Ky 
Kı 


Increasing the magnitude of C, will decrease both к. 
апа 1/K, resulting in the ce; cross-coupling 
effects on the angle-of-attack (a) and the influence of °D 
in the a. 

In the section III the equation of NE derived: 

LETS = (-as |w) i CRIE N) 6:5) 
which involve (a). Therefore, the reduction of inertial 
and kinematic effects on (a), results in the reduction of 
the effects in the achieved body-fixed acceleration EU 

In the nonlinear studies the equation for the achieved 
inertial acceleration was derived as: 

Му = "ag - co5 + Ny, -SIN (V.C.6) 

It is apparent from the equations (V.C.4), (V.C.5) and 
(V.C.6), that an increasing of Са results in the reduction 

a 


of kinematic and inertial effects (i.e transients in the 


response) in the achieved inertial acceleration (пу), which 
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it is desired in order to meet the requirement of time 
constants T3) and percents of overshoots and 


undershoot, as referred in Appendix C. 


D. MINIMIZATION OF THE KINEMATIC AND INERTIAL 

CROSS-COUPLING EFFECTS IN ELLIPTICAL AIRFRAME 

As discussed in the nonlinear studies, using the pitch 
angular rate q as feedback and under the influence of the 
inertial and kinematic effects, the achieved inertial 
acceleration (n3) for both airframes doesn't meet the 
requirements of the time constants, overshoot and 
undershoot, denoted for the linear autopilots, referred ta 
in Appendix C. 

Figure 4.5 and Table II show that the achieved 
inertial acceleration (nj) For elliptical airframe with 
feedback the pitch angular rate (q) has overshoot 11.25% 
and в = 0.48 sec, iS 0.41 sec. 

Figures 4.48 and Table III illustrate that the 
achieved inertial acceleration (nz) for circular airframe, 
with feedback the (q) has overshoot 5% and T=-,46 Sec, 
бо = -.62 sec. Непсе the (n,) £or both airframes is not 
acceptable. 

In order to minimize the inertial and kinematic 
cross-coupling effect, namely to have smoother transient 


response of n; which will meet the requirements of time 
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constants and overshoot, the magnitude of Са was increased, 
a 


utilizing as feedback the angle-of-attack and its 
derivative. 


In order to increase the magnitude of Ст the 
a 


stability augmentation theory was used. 


Considered as stability derivative augmentation = 


and > increase the magnitude of M, and M; directly, as 
referred in [Ref. 6], the increases in the Ma is 
7 


essentially an increase in Сп and hence the inertial and 


a 
kinematic effects on the angle-of-attack are reduced. 


In the pitch control law (Figure 3.1) a SAS exists 
using as feedback the pitch angular rate (q). In order to 
design the SAS with the angle of attack (a) as feedback, 
the root locus and lead-lag compensation techniques were 


used. 


1. Stability Augmentation System (SAS) With Feedback 
of the Pitch Angular Rate (q). 
In Figure 5.1, the stability augmentation system 
is shown, with feedback of pitch angular rate q where the 


transfer function of the q/sS, for the uncoupled pitch 


Р 
dynamic model, given in section II.C. The closed-loop 


transfer function of the above SAS is: 
K (=+1) (5.3) (-41.35- 5 - 7.7118) 


5 (mg 1) (40455-6239) 
а5--------------- ШЕБІ! 


(toma (+4) (53-3) (-41.35- 5-7-7418) 


s (513) (5 +1) (5+0. 0815-11299) 






әшелултү теотадатття 3jO LLgO '[euueuo чозта jo 


а 
9/b uyIM (SWS) waysAs чотзезиэшБиуч Азтттаез5 





Bap) b | 
(Sap) GETIT -S T810 tS 
IFA LS ISE TH 
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Manipulating equation V.D.1-1, it results: 
E х (asegar s- 242.563 -$ ~ 444-8962 ) 

0.0053. 5'+1.0004-5++0.11-5- 1128-5 + К(-5.1695- 42.3245 - 3-3418) Сол, 
The poles and zeros of the SAS are as follows: 
ив. = 183.43 
b. P) = -3.45 


Es P3 = 3.269 


ar РА = 0.0 
апа 

е. 21 = -—8.0 
f. 2. - -(.186 


The other two zeros are at infinity. 

The root locus of the system is shown in Figure 
5.2). 

The K was selecting at the value of —3.07, because 
at this value, the system has two dominant roots at Piro = 
-9.212%;7.606 which dominate the response of the system, 


having the following characteristics: 


Y -{ 7.606 
Bang te en 


b. 3 = cos Ə = 0.133 


ce (j= Us V 4-35 = ži 1-32 = 7645 («од ес) 


d. Uh = = = 11.203 (vadIsec) 


= 427$ 





where: 
с : damping constant 
ш $ frequency of system 


ш: natural undamped frequency 
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2. SAS with Feedback of the Angle of Attack (a) and 


mate OL angle of attack (a) 
The SAS with feedback the pitch angular rate q 


(Figure 5.1) and K = -3.97, is stable and the dominant 
roots gives damping ratio < = 9.73 in the response of the 
system. 

To minimize the coupling effects the angle of 
attack (a) was selected as feedback in the previous SAS. 
It is desired the SAS to have the same behavior as before, 
for purposes of stability and response of whole the pitch 
channel. 

Pole zero cancellation technique is used, to 
obtain the same poles and zeros, as in the case with 
feedback the pitch angular rate q. 

For above purposes, a compensator was used in the 
SAS, aS shown in Figure 5.3. 


The transfer function of the compensator is: 


_ -41.354-5 - 4.1418 
G(s) = -0.02.5 - 41.35} 


The closed loop transfer function of the SAS 15 as 


Wen) 


follows: 
sn). PEREA c 34.5 \- сс аен ) 
БЕ о аә LEE ^^" "ершк ee P were 5 
Mis) - 0.02-$- 01.253 $ apri 5°+ 0.1845 - 11.279 


44 B Ur =) 53.3 (ашын ) 
545 |) -0.02-5-4И.35% 5 at) 840.181-5- 11.273 
(V.D.2-2) 


where the transfer function a _ -0.02.$-H1.351 is 


бр $4+0,484.5-11119 


derived in the section II.C.3. 
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Rearranging and manipulating the equation 


(V.D.2-2), it becomes: 


6(5)= 


апа 


dominant 


3. 


2 
к (5.9234 564229,9 5+ 400292.449:5 + 13235.09) 
-0.0001-55% 0,2215 54.41.399 5° 4.39 52. 66.4665 + К (0.1034 5+ 214.6 54-1150.15 6319.9) 


(V.D. 2-3) 


The poles and zeros of the system are: 


ЕЕ, “1877.931 
P, = -207.391 
P4 = -3.447 

Pa = -3.271 

P. = 0.0 


2) - -2067.838 

2 = -8.0 

2: = -0.186 

The other two zeros are at infinity. 


The root locus of the system is shown in figure 


Selecting the K = -3.97, the system has the same 
roots P} 5 = -8.249257.577, as in section V.D.1. 
г 


Acceleration Loop 


Considering the achieved acceleration loop in the 


Pitch channel, the modification of the SAS, effects the 


transmission gain of the above loop. 


коз 
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In order to compensate the effect in the 
transmission path, a compensator was added as shown in 
Figure 5.5. 

The transfer function of the compensator is the 
inverse of the V.D.2-1: 

1 -0.02- S- 41.357 
EISE eT (V.D.3-1) 
бс)  —-44.351- S - 7.4113 
Modifying in this way, the transmission path of 


the acceleration loop, the gain of path, is the same as in 


Ehe case of par. V.D.l. 


4. Measuring of the Minimization in the Kinematic and 
Inertial Cross Coupling Effect 


In order to measure the minimization a CSMP 
program (Appendix J) was written using the equations of the 
nonlinear aerodynamic models, the equations of nonlinear 
lateral control law and the equations of the modified pitch 
control law. 

The equations of the modified pitch control law 
are the following: 

a. Acceleration lag network 

Х= -150х+ 150 N?B 65520) 

Ы. Acceleration second order compensator (Figure 
5.5) Y _ €0.008). C0.02- S -ut.353) 

(1110, -X) (£ +1) (— 44.359-5-7.7118) (V.D.4-1) 
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Utilizing inverse Laplace transformation and 
rearranging the above equation, one has: 
Y 46.1961 Y +483 У = -0.000232 (5.11 0, -X)- 0.«186- (1-149, - X) (V.D.4-2) 
Using space representation of a system, in which 
the forcing function involves derivative 


[Ref. 5, pp 675-678], one obtains: 


№ = W,-0.00025751, +0.000232 - X (V.D. 4-3) 
Ño 2 -1.4881 W - 6.4867-W, -0.5312-M,HOMFE6X (V.D.4-4) 
eom (V.D.4-5) 


с. Compensator of Actuator 
See (+4) 41.3515 -1.1118) 
(М- 24533) | 5 (-0.025-41.35%) ашы 
Using the method for state representation and 


rearranging the equation УСЫ 4—5 yields: 


бр = Sp +K: 532.382 x10 (Y-al5t.3) (V.D. 4-6) 
OP, = -206}.35 др» - К. 1100. 836 x40 (y-al57.3) (V.D.4-7) 
Ôp = de - k- 258.48 (Y-als1.3) (V.D.4-8) 


а. Equation of Actuator 
бь- -18%4.9р-і. 4095,3 -бр, (ІУ.В.1-4) 
Selecting K23.07, for the compensator of SAS, with 
feedback (a) the time constant i! (i.e time constant of 
achieved maneuver acceleration due to the first command) is 
9.45 sec, the time constant т, (for second command) is 8.41 


2. 


Sec and the undershoot is 20.73 percent as shown in figure 


5.6. 


203 






ES lg тать 
Teot23d1113 JO 1199 “Тәчиецо цоҙта perjtpow 


(3) eut; SA (Zu) чотаеләтәооу тетіләиі peAetuov 9°с “Бід 


(22493411 
0078 0274 оћ'9 09'5 08*h 00*h 0276 оћ'а 0971 0870 00° 
пио А љиљани] аи 1 





09 eS 


oh°e- 00 °Е- 


0274- 0871- 
(SÓ ZN 


09°0- 


00°0 


гы 
ж 
0880 


043931 


og'i 


204 






The above response of the inertial achieved 
maneuver acceleration does not have, smooth transient, 
because of the large undershoot. 

Due to the above reason, another K was selected 
equal to -19.37, for minimization of the coupling effects. 

With this selected K the system has the following 


poles, shown in figure 5.4 


Р} ‚2 = -98 596.239 
Pz = -1879.06 


with the same zeros as section V.D.2, the 
dominant poles of the system are P,,9, and the 


Characteristics of the response are: 


а. Ө = 44.24 
DR g = 0.716 
С. w = 93.775 rad/sec 
а. ш = 134.412 rad/sec 


It is shown in known (Ref. 6) that the ш сап бе 


expressed as: 


Ma 
1 





Wy 


1 


(У.р. 4-9) 
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and 


Maz Cm ŽL. PMA RUNDE 0 


where: 

Y: specific weight of air 1.4 

I: moment of inertia 

P: Local pressure 

M: Mach number 

A: Area 

Comparing the O, = 134.412 for the K = -19.37, it 
is apparent that it is much larger than the D 11.203 for 
К - -3.07 and hence according to the equations V.D.4-9 and 
V.D.4-10, the C. increases essentially, resulting in 
minimization of atis and inertial cross coupling 
effects. 

Figure 5.7 shows, that with K = -19.37 and 
feedback in SAS the (a), the achieved inertial maneuver 
acceleration has Um -.49 sec, т, = -.38 sec and overshoot 
about 0%. 

The above analysis illustrates that the selection 
of K = -19.37, instead of K = -3.07, keeps the damping 
ratio about constant ¢ = 0.71 (Сс = 0.73 for K =-3.97) and 


increases the Wr which essentially increases the C. ; 
а 


resulting in very smooth response of the achieved maneuver 


acceleration. 
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Using the same particular compensators in the 
Pitch channel shown in Figure 5.5 and feedback the rate of 
angle of attack (a), the Pitch channel becomes 
uncontrolable for both selection of K = -3.097 апа К = 
-19.37, because the pitch tail incidence exceeds the limits 
оғ 319°. 

Also, selecting the K = -19.37 in SAS of Figure 
5.1 with feedback the pitch angular rate а the Pitch 
channel becomes uncontrolable again, due to the exceeding 
Sethe limits іп $. 

Therefore, for minimization of the kinematic and 
inertial effects (i.e improving of the transient response 
of the achieved maneuver acceleration) for this particular 
airframe and flight conditions, the C, was increased, 
using the aforementioned compensators 2 ehe piten control 
law (Figure 5.5) and feedback the angle-of-attack (a). 

Comparing the achieved inertial acceleration (п2) 
with feedback (q) in the Pitch channel (Figure 4.41) and 
the achieved inertial acceleration (nz) with feedback (a) 
and K= -19.37 (Figure 5.7), it is concluded that for these 
particular airframes and flight conditions the feedback of 
angle-of-attack (a) is more desirable, than the feedback of 
(а) and (a), because in the achieved inertial acceleration 
has been reduced the kinematic and inertial effects, 


meeting with this way the requirement of the time constants 


and overshoot, denoted in the linear studies. 
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E. MINIMIZATION OF KINEMATIC AND INERTIAL COUPLING IN 

ELLIPTICAL AIRFRAME 

Making a comparison in Figure 5.8, which shows the 
kinematic cross-coupling (-p8) into а, in the pitch channel 
with feedback pitch angular rate (g) and Figure 5.19 which 
shows the -p8 in the modified pitch channel (Figure 5-5) 
with feedback the angle-of-attack, it is concluded (Table 
IV) that the kinematic coupling is increased from 
(-09.988-59.877) to (-29.124-79.383), when the kinematic 
coupling effects are reduced by the modified pitch channel. 

Figure 5.9 shows the inertial coupling (pr) into q in 
the nonlinear Pitch channel, with feedback the pitch 
angular rate (q). 

Figure 5.11 shows the inertial coupling for the 
modified Pitch autopilot with feedback the angle of attack 
(a) (Figure 5.5). 

Comparing the above figures, it is concluded (Table 
IV) that the inertial coupling is increased from 
(-1°.692-168°.317) to (-29°.373 - 181°.415), when the 
inertial coupling effects are reduced by the modified pitch 
channel discussed in previous section. 

Hence, the minimization of the kinematic and inertial 
effect does not mean necessarily minimization of the 


cross-couplings themselves. 
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For minimization of the kinematic and inertial cross 
coupling, the same modified nonlinear pitch control law was 
used as that of section IV.D (Figure 5.5), except the 
acceleration compensator which was removed from the 
transmission path of the n. as shown in Figure 5.12. 

For purposes of analysis a CSMP program (Appendix J) 
was written, using the following equation for the nonlinear 
pitch control law and the same equation for lateral control 


law and dynamic models, as in section IV.B and IV.C. 


= — 150Х + 150. 1:6 (IV.B.1-1) 
Y = -6-Y- 0.5305 Nc + 0.43-X (IV.B.1-2) 
др, E ӛр, + к. 532. 382 X 407. ( Y- ad51.3) (V.D.4-6) 
dp, 2 -2063.95 -OP> – К: 1400.956 x10° (y-als}.3) (V.D.4-7) 
др, = А ~ к. 258.43 (Y- alst.3) (V.D.4-8) 
dp = -188.4. Ор + 10395.3 · дрс (ІУ.В.1-4) 


where K = 19.37. 

Figure 5.13 shows the minimized kinematic coupling and 
Figure 5.14 shows the minimized inertial coupling. Figures 
5.13, 5.14 and table IV illustrate that the minimization of 
the kinematic and inertial coupling is 264% and 195.26% 
respectively from its values in the nonlinear CBTT 
autopilot for elliptical airframe. Figure 5.15 shows the 
achieved inertial acceleration (no), with the Pitch channel 
modified for minimization of kinematic and inertial 


coupling (Figure 5.12). From this figure it is concluded 
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TABLE IV 


KINEMATIC AND INERTIAL CROSS-COUPLING FOR 
ELLIPTICAL AIRFRAME 


CONDITION INERTIAL (deg/sec”) KINEMATIC (deg/sec?) 


MIN MAX MIN MAX 
Elliptical airframe 


with cross- 
coupling -1.692 168.847 -0.9881 52877 
(section IV.C) 


Elliptical ideal 
airframe (section -1.698 167.37 -2.124 7.383 
IV.E) 


Elliptical airframe 

with modified Pitch 

autopilot -20.373 181,415 -1.88 8.644 
K = -19.37 

(section V.]) 


Modified Pitch 
autopilot K = -30.4 -7.263 169.371 -1.956 8:71 
(section V.D.) 


Modified Pitch 

autopilot for -0.0577 57.5322 - 0.523 1,363 
minimization of 

cross couplings 

(section V.E) 

K = -19.37 


Note: K is the gain of the added compensator in SAS of the 
Pitch autopilot. 
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that the kinematic and inertial effects (i.e transient in 


achieved acceleration response) is increased. 


Е. MINIMIZATION OF THE KINEMATIC AND INERTIAL 

CROSS-COUPLING EFFECTS IN CIRCULAR AIRFRAME 

As shown in Figures 4.28, 4.40 and 4.42 the achieved 
maneuver acceleration in the circular airframe has very 
smooth response, in regards to the overshoot or undershoot 
en it. 

The effect of the inertial and kinematic 
cross-coupling in the nonlinear aerodynamic model is that 
the response of the inertial achieved acceleration has a 
slowing transient starts at 3.5 seconds, as discussed in 
section IV.G. Непсе, minimization of the effects means in 
the circular airframe improvement of time constant during 


the second guidance command. 


l. Stability Augmentation System (SAS) with Feedback 
the а. 


Using the same method for minimization, as in 
section V.D, the SAS of the pitch channel for the circular 
airframe with feedback q is, as shown in Figure 5.16. The 
transfer function of the q/ Š, is given in section II.C. 

The closed-loop transfer function of the above SAS 


15: 
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(5%) 53.3 (- 59.103885 – 6.9845) 
(+4) (1+4) (5% 0444391 5+ 44.33 ) 
26) = ИИС л Оуу (54 04484 4453) — Ener (отима) ок. У) 
к. [255+1) (5+3) (.53.1083-5 - 6.9845) 


63 (e +) (a +1) ($+0.1478.5 +44,3346) 


Rearranging equation V.F.l-1, one obtains: 


. K (20.38-5^4-30200.25 $4 203826 S t 23612.9 ) 
e IT 
0.037-5"+ 7.004 $° + 3.679-S° + 310.39 :5 +24.33 +к. (-9.024:5- 60.199-3- 6.479) 
(V.F.1-2) 


2(5) 


The poles and zeros of the SAS are: (1) 
P,= -188.438; (2) P, 3 = НВО порт В, = 0.143; 
(4) ЕЕЕ -6.55 апа (5) 22 = -0.118. The other two zeros 
are at infinity. The root locus of the SAS is shown in 
Figure 5.17. The selected value of K is -15.6 and the 
dominant roots are Pi 2 = —19.7361#3j)7.536 (Figure 5.17). 
With those dominant roots, the response of the system has 9 
equals to 35.06 degrees; & equals to 9.818; w equals to 


7.545 rad/sec and 9. equals to 13.125 rad/sec. 


n 
2. SAS with Feedback of the a 
For minimization of inertial and kinematic effect 
in the circular airframe, the same technique was used (i.e 
pole-zero cancellation technique) as in section V.D.2 and 
the resulted pitch control law with feedback of 
angle-of-attack (a), is shown in Figure 5.18. 


The transfer function of the added compensator in 


the SAS is given by equation V.D.2-2. 


222 (V.D.2-2) 
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The closed loop transfer function of the referred 


before system is given by equation V.D.2-3: 


5 
(ся?) .[2591:5-6:93 \ ү 51.3 X) Ea ) 
-00395-54 / |Жаз!) Қ 5%0448544,33 | 












Ko (+4) 513- (-5945- 6.99) (-0.039 5- 59.1) 

a 54.3 (Spy +4)- -0.039-5-59.1) (S741) (s4o.t48-5 444.33) 
(V.D.2-3) 

where the transfer function ars, is given in the 

section II.C.3, rearranging and minor manipulating the 

above equation, one obtains equation V.D.2-4: 


209- 038 3+ 30400.25 5+ 203825.65+23642.3 ) ET 
(V.D.2-4) 
The poles and zeros of the modified SAS are, 
Pj ^ -1499.623, P, - -188.53, ый 2530-0560, DES 
-0.142, 2) = -1499.6, 25 = 6.55, Z, = -9.1182 and all 
other zeros at infinity. 
The root locus of the system is shown in figure 
5.19, it is apparent from above Figure, that for k=-15.6 
the dominant roots are the same as those in section V.F.1. 
Considering the transmission of acceleration nc 
path, as in the elliptical airframe, a compensator was 
added in the path, for compensation of the effect in it, 


which is created by the added compensator in the SAS, as 


shown in Figure 5.28. 


225 





-220.00 


Fig. 





REAL AXIS(UNITS PER INCH! = 20.0000 
IMAG AXISCUNITS PER INCH! = 19.0000 


BOGEOCZBIPOPBNRHETSTIS- BF EBT? AUTOPILOT 
CIRCULAR AIRFRAME 








20.00 -100.00 
REAL AXIS 





=200.00 -190.00 -180.00 -140.00 ~t -80.00 


+ 


5.19 Root locus of SAS with a/3$ 
CBTT of Circular airframe 


226 


р 









4D. 00 
IMAG AXIS 


.00 





эшелулте летподто 103 Тэцчиецо цозта рэтутрои 0276 “Бтд 











80}'bS - $ h6£0'0- 
h8b'9 - $80b'bS- 


+ h8b'9 - $gorbs- 














Ct hh $3EhJ'0-s 
80465 - $ hbeo-o- 


227 





В = Measuring of the Minimization in the Kinematic and 
Inertial Cross-Coupling Effect 


A CSMP program (Appendix K) was written using 
equations of nonlinear lateral control law (section 
IV.G.2), of the aerodynamic models (section IV.4) and the 
following equations of the modified nonlinear pitch control 
law: 

a. Acceleration lag Compensator 

X 2 -150X + 150 Үр (11.64) 

b. Acceleration Second-Order Compensator (Figure 


3.20) с 
B: m K: (5-4) (-0.03945-594) «ye 3-1) 


(ч- X) 5 (5+1) (- 59.103. 5 - 6.9845) 
Utilizing inverse Laplace transformation and 


rearranging the above equation, yields 
V, 15.1133 Y 0,591 Y z-4.0.02283 (5, - X) - 3.81 (8 -X)- S (nex) 


UE 3-2) 
Using state-space representation referred to in 


(Ref. 5: pp 675-678), one obtains: 


W, = Wa - К, - 0.02288 (ме-Х) (V.F.3-3) 
W, = Ws - Ka - 34.850 (n, -X) (V.F.3-4) 
W; = -0.59412. W. - 5.118. + Ka- 473,39 (n,-X) (V.F.3-5) 
Үз = Wy (V.F.3-6) 
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с. Compensator of Actuator in SAS 
де а (+) (-59.1088.5-6:989) (v.F.3-7) 

(Y3- |543) Ё (= eH) (-0. 0334-5- 59.104 ) 

Utilizing, the aforementioned method for state 
representation and rearranging the equation V.F.3-7, it 
yields: 

бр, = др, + к (43806.0838). (Y3-alst3) (У.Е.3-8) 

Р› =-249.433 вр - 1500 др,- х.4.3 40 (равз)(У.Е.3-9) 

дрес = др, _ Ko (32.744) (Уз- а[53.3) (V.F.3-19) 

d. Equation of Autuator 


(IV.G.3-8) 


др = -188.4 Әр 10495,5. bp 


Selecting the K = -15.6 and Ra -0.0387 
(Figure 4.26) and feedback the pitch angular rate (q), the 
achieved inertial acceleration n,, has Tu = -,34 sec, т. = 
0.62 sec and overshoot 5$ (Figure 4.40, Table III). 


This response of the n, is unacceptable 


2 
because it doesn't meet the requirement on the time 
constant, referred to in Appendix C. 

Selecting K = -15.6 and Ka = -0.0387 (figure 
5.20) and feedback the angle of attack (a) in the pitch 
channel the achieved inertial acceleration n; (Figure 5.21) 


has E - -.35 sec, т = @.6 sec and overshoot 7%, which is 


again unacceptable. 
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It iS known (Ref. 6) that: 


an = EA 
Wy Y my ЕРМ?.А 
Hence, in order to reduce the т, it is needed 
to increase the Cn ° Selecting K = -70.88 and Кл = -0.0337 
(Figure 5.29) aan о. the angle of attack (a), the 


(MS) 


dominant roots of the SAS (Figure 5.13) are Pira = 
=99.553j87.95 and the characteristics of the system 
response are 9 = 44.16 degrees, ¢ = 0.79 , w = 87.932 


126.212 rad/sec. With that selected k the 


rad/sec and e 


Z decreases from 0.818 (when K= -15.6) to 0.72 апа the ЈЕ 


increases from 13.125 rad/sec (when K = -15.6) to 126.212 
rad/ sec. 
Figure 5.22 shows, the achieved inertial 


acceleration n when K = -70.88, K, = -9.0387 and feedback 


, 
2 A, 


the (a). The e = 9.33 sec, е 0.48 sec and overshoot 


14.33 %. It 1S apparent that the ES decreases, which is 
desirable but the overshoot increases, because the damping 
ratio ¢ decreases. 

In order to reduce the overshoot and to keep 
the time constants below the 0.5 sec, the Кл reduces to 


-0.0274 (Figure 5.29) and the gain K (Figure 4.27) of the 


YE 


coordination branch in the nonlinear lateral control law 
increases from 9.458 to 1.8. 


Using K = 1.9 (Figure 4.27) in the 


YE 


coordination branch in the lateral control law, K, 
1 
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= -0.0274 in the acceleration compensator and K = -79.88 in 
the added compensator in the SAS with feedback the 
angle-of-attack (a) (Figure 5.20), the obtained achieved 
maneuver acceleration has ul = -,43 sec, т. = = 37 sec 
overshoot 8.9%, as shown in figure 5.25 and table У. 

Using in the pitch channel with the 
aforementioned particular compensators, as feedbback the 
rate of angle of attack (a) and making all the above 
K and K 


selection of K the whole system becomes 


үр 


uncontrolable because the required tail incidence angle 


А’ 


exceeds the limits of 2109, 


Also using the selections of K, K K and 


А. "92р 
feedback pitch angular rate q (Figures 4.26, 4.27), again 
the whole system becomes uncontrolable, because the 
required tail incidence angle exceeds the limit of Eng. 
The above analysis illustrates that the 
angle-of-attack is more desirable feedback than the (q) and 
(a), in the nonlinear pitch control law for these 
particular circular airframe and flight conditions, because 
it reduces the kinematic and inertial coupling effect, such 
that the achieved inertial acceleration (nz) meets the 
requirements of the time constants and overshoot, referred 


to in the linear studies. 
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TABLE V 


CHARACTERISTIC OF THE INERTIAL ACHIEVED ACCELERATION (nz) 


FOR CIRCULAR AIRFRAME 


CONDITION 


K E 0.0387, K = -15.6 


-15.6 


ж 
| 


К, = 0.0274, 


Е 09-0274, К = -70.88 


Ty (sec) T (sec) overshoot 
0. 35 0.6 7 

ЖЕК 0.48 14.33 
0.435 DS 8.6 

0.43 o. 37 8.9 


Note: K,: pitch autopilot acceleration gain in acceleration 


compensator 


ту? 63 percent time constant of achieved maneuver 
plane acceleration due to first command 

T: 63 percent time constant of achieved maneuver 
plane acceleration due to second command 


autopilot 


ур: gain of the coordination branch 
К: gain of the added compensator is the SAS of pitch 
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G. MINIMIZATION OF KINEMATIC AND INERTIAL COUPLING IN 

CIRCULAR AIRFRAME 

As in the elliptical airframe comparing the values of 
the inertial and kinematic cross-coupling for the circular 
airframe, namely, the (P.r) and (-p8), referred to in table 
VI, it is concluded for this particular analysis that as 
the coupling effects in the nonlinear Pitch autopilot are 
reduced, the cross-couplings themselves are increased. 

In order to minimize the cross-couplings, the pitch 
control law was modified, as shown in Figure 5.26. 

For purposes of analysis, a CSMP program (Appendix K) 
was written, using the following equations for the 
nonlinear modified pitch control law, lateral control law 


and dynamic models, as in sections IV.G and IV.H. 


o 


2X = -150X +150 -Nig oen 
2. W= Wa - X4: 0.02288 (XY (V.F. 3-3) 
3. Wo = W3 - Кл. 34.850 (Nc-X) (V.F.3-4) 
4. Wa = -0.594 \№ - 5.4184 №5 +КА: 13.385 (ух) (У-Ғ.4-5) 
5. Ya= Wi (V.F.4-6) 
є. 85 =-0443.6р -0.34058 (5-@.13)-22308 -(%-0453.3) 
(V.E.4-7) 


which is the equation IV.G.1-4, except of the 
feedback, which is the angle of attack (a) instead of the 
pitch angular rate (q). 


7. dp = -4 33.4 - dp - 10395,3-4рс (IV.G.1-5) 
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TABLE VI 


KINEMATIC AND INERTIAL COUPLING FOR CIRCULAR AIRFRAME 


CONDTION 


Circular airframe 
with cross- 
couplings 
(section IV.G) 


with cross- 
coupling and 
faster response 
K. = -0.0387, 

К = -15.6 
(section V.F) 


with cross- 
coupling and 
faster response 
БА = -.0387, 

K = -70.88 
(section V.F) 


with cross- 
coupling and 
slower response 
Ka = -0.0274, 

K = -15.6 


(section V.F) 


with cross- 
coupling and 
slower response 
Ka - -0.0274, 

К = -70.88 


(section V.F) 


with cross- 
coupling 

slower response 
and increased gain 
in the coordinated 


branch K = 1.8, 
К = -70.88 


(section V.F) 


with cross- 
coupling slower 
response and 
modified pitch 
control law as 
section V.G 


Note: 


INERTIAL (deg/ sec?) 


MIN 


-0.008 


-0. 0768 


-0.90212 


= 050107 


-0.0743 


-6.0791 


-2.264 


МАХ 


92.794 


199.573 


1932533 


118-294 


109.411 


199.841 


38.386 


239 


KINEMATIC (deg/sec^) 


MIN MAX 
INNO 2.012 
-2.04 0.0355 
Els 0.0436 
=2- 97 2.0387 
=2 32 9.0431 
-0.8662 1.0541 
-0.00788 -0.6021 


КА: К, Кур as referred in section V.F and V.G. 
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Figures 5.27 and 5.28 show the minimized inertial and 
kinematic cross-coupling, with feedback the 
angle-of-attack, with КА, -0.0274 and К - -70.88. 
Comparing the values of the above minimized inertial and 
kinematic coupling, with the values of the Kinematic and 
inertial coupling when the pitch channel has feedback the 
pitch angular rate (g) and Ka = -0.0274, as given in table 
VI, lt is concluded that the couplings when the pitch 
channel has feedback angle-of-attack (a) and the 
aforementioned particular compensator (Fgure 5.26), were 
reduced 228.293 for inertial and 19093 for kinematic. 

Also, comparing the Figure 5.25 which illustrates the 
achieved inertial acceleration n, for Circular airframe 
with minimized coupling effects and Figure 5.29 illustrated 
п, for circular airframe with minimized coupling, it is 
concluded that as the coupling reduces, the coupling 


effects increases. 


H. CONCLUSIONS AND RECOMMENDATIONS 

The conclusions that follaw are based on a single 
representative flight condition (Mach 3.95 and 60 kft 
altitude) and for the particular analysis which was done in 
this chapter. 

1. From the studies in this chapter, it is concluded: 


a. Transients in maneuver plane acceleration are 
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caused by kinematic and inertial coupling between pitch and 
yaw dynamics through missile roll rate. Transients are in 
the form of overshoots and undershoots (section V.B). 

b. Using as feedback in the nonlinear pitch 
channel for both airframes the pitch angular rate (q), the 
transients in the maneuver plane exceeds the denoted 
requirements and the achieved maneuver acceleration is not 
acceptable (sections V.D and V.F). 

c. Using as feedback in the nonlinear pitch 
channel for both airframe, the angle-of-attack (a) and 
modifying suitably the pitch control law (Figure 5.5 and 
Figure 5.28), the transients are minimized below the 
desired requirement and the achieved inertial acceleration 
(nz) is acceptable. (sections V.D and V.F). 

а. Using as feedback the rate of angle of attack 
(a), the whole system becomes uncontrollable because the 
pitch tail incidence ($5) exceeds the limits. (sections 
ПО МЕ). 

e. Using as feedback the rate of angle of attack 
(a) and modifying suitably the pitch control law (Figure 
5.26), the kinematic and inertial coupling are reduced 
dramatically, comparing with their values when the used 
feedback is the pitch angular rate. (sections V.D, and V.G) 

f. Using as feedback the rate of angle of attack 


(a), in the modified pitch channel (Figure 5.26) for 
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minimization of kinematic and inertial coupling, the whole 
system becomes uncontrolable, because the pitch tail 
incidence ( $5) exceeds the limits. (sections V.E and V.G) 

g. Minimizing the inertial and kinematic 
cross-coupling effects, the values of kinematic (-p8) and 
inertial (pr) coupling increase. (section V.E and V.G) 

h. The above conclusions show, that for the 
particular airframe, flight conditions and analysis is not 
possible to achieve minimization of the kinematic and 
inertial cross-coupling effects and simultaneously 
minimization of the cross-couplings themselves. 

2. Recommendations 

Considering the studies of this work and the 
conclusions, it is recommended: 

a. Further analysis of the minimization of the 
kinematic and inertial cross-coupling effects and the 
cross-coupling themselves in other airframes, and with 
other analysis to achieve simultaneous minimization. 

b. Further analysis of the CBTT autopilots in 


other flight conditions for verification of the above 


results. 
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APPENDIX A 


MISSILE SIZING AND MASS PROPERTIES 
In order to provide a realistic missile based on the 
configuration concepts tested aerodynamically in [Ref. 1], 


the models have the following characteristics. 


Characteristic Circular Elliptical 
Length (in) 168 168 
Max. Diameter 24 

Max. Major axis (in) 414,57 
Max, Minor axis (in) 13:86 
Weight (lbs) 2525 2475 
I,, (slug-£t“) 40 110 
b. (slug-£t) 804 790 
I, (slug-ft^) 810 853 
c.g distance from L.E (in) 100.8 100.8 
Reference length d (ft) 2 2 
Reference area S(ft2) > а 
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APPENDIX B 


NONLINEAR AERODYNAMIC DATA 


The aerodynamic data were taken from [Ref. 2]. For 
reference, the normal and pitching moment curves at M = 
3.95 have been produced in fig. B.1 and B.2 for the two 
configurations in [Ref. 2]. The aerodynamic derivativies 


Су’ Ch! Cor with respect to sideslip angle 8, yaw control 


and roll control бы аге presented in figures B.3, B.4 


бу» 
апа B.5 of [Ref. 2] as used in the computer simulation, 
namely, as piece wise linear segments for ease in 


interpolation. 
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APPENDIX C 


REQUIREMENTS FOR UNCOUPLED AUTOPILOTS 


The requirements for the classical design technique of 
the uncoupled autopilots are given in [Ref. 2] and are the 
following: 

1. Acceleration time response 
a. 63 percent time constant of 9.5 seconds for a 
step command of acceleration at the flight condition of 
interest, (i.e M = 3.95 and altitude 60kft) and small 
angles-of-attack. This response is representative of a 
tactical missile of this size. 
b. Overshoot 10 percent 
с. Zero steady state error in acceleration to 
reduce variations of quidance navigation qain. 
2. Relative stability 
Gain margins 6db, phase margins 30 deq with a 
goal of 12 db and 50 deg. 


3. High Frequency Attenuation in Actuator Command 


Branch 

It must be 5, 15 db at 100 rad/sec and zero 
angle-of-attack. This requirement will provide sufficient 
high frequency attenuation for 30 Hz actuator and for body 


bending modes, when high frequency filters are added. This 


requirement also limits autopilot speeds of response. 
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APPENCIX D 
FOR PITCH LNCOUPLED 
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APPENDIX E 
СМР PROGRAMS FOR YAw UNCOUPLED AUTOPILOT 
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APPENDIX F 
CSMP PROGRAMS FOR RCLL UNCOUPLED AUTOPILOT 
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APPENDIX G 
CSMP PRECRAMS FÜR LINEAR CBTT AUTOPILCTS 
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